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Abstract 

We propose a theoretical framework for R-parity violation. It is realized by a class 
of Calabi-Yau compactification of Heterotic string theory. Trilinear R-parity violation in 
superpotential is either absent or negligibly small without an unbroken symmetry, due 
to a selection rule based on charge counting of a spontaneously broken U(l) symmetry. 
Although such a selection rule cannot be applied in general to non-renormalizable opera- 
tors in the low-energy effective superpotential, it is valid for terms trilinear in low-energy 
degrees of freedom, and hence can be used as a solution to the dimension-4 proton de- 
cay problem in the minimal supersymmetric standard model. Bilinear R-parity violation 
is generated, but there are good reasons why they are small enough to satisfy its upper 
bounds from neutrino mass and washout of baryon/lepton asymmetry. All R-parity violat- 
ing dimension-5 operators can be generated. In this theoretical framework, nucleons can 
decay through squark-exchange diagrams combining dimension-5 and bilinear R-parity 
violating operators. B — L breaking neutron decay is predicted. 
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1 Introduction 

In supersymmetric extensions of the standard model, renormalizable operators 

W 3 XLEL + X'LQD + X" DUD 5105 (1.1) 

break baryon and lepton number symmetries, and hence a proton decays too rapidly. Either 
the coefficient A" of the baryon number violating operator DUD or A and A' of the lepton 
number violating LEL and LQD have to be highly suppressed. The most popular solution to 
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this dimension-4 proton decay problem is to assume an unbroken R parity (or matter parity), 
which removes the operators in ( ll.ip altogether. 

Various alternative solutions have also been discussed in the literature, some of which are 
found in a review article [1]. Some solutions assume discrete symmetries other than R-parity 
(see also [2] and references therein), so that either the first two operators or the last one 
are (is) forbidden by the discrete symmetry. Although phenomenological consequences of these 
solutions are quite different from those of R-parity preserving ones, these two classes of solutions 
share one thing in common, an unbroken discrete symmetry. In Calabi-Yau compactification 
of Heterotic x Eg string theory, however, discrete symmetries are found only at special 
points (or subsets) in moduli space. Reasons are not clear why such vacua have to be chosen. 
References [3lll], for example, clearly expressed dissatisfaction to solutions relying on unbroken 
discrete symmetries. 

This article presents an alternative solution to the problem without assuming an unbroken 
symmetry. The essence of the solution is an extra U(l) gauge symmetry with a Fayet-Iliopoulos 
parameter. The U(l) symmetry is broken spontaneously at high energyjl] allowing for large 
Majorana masses of right-handed neutrinos. No unbroken symmetry is left at low energy, but 
its legacy still remains. There is a selection rule [5] (also known as SUSY-zero mechanism) in 
how the U(l)-breaking vacuum expectation value (vev) can appear in superpotential of low- 
energy effective theory, and this rule may be used to make sure that the trilinear R-parity 
violating couplings (11.11) are absent. In this sense, the solution in this article is certainly along 
the line of models in [6l[7l[8l[9]. It should be reminded, however, that there are many subtleties 
in how to use the selection rule in low-energy effective superpotential. It is one of the purposes 
of section [2] of this article to clarify when the selection rule can be used and when it cannot. This 
solution fits very well with Calabi-Yau compactification of the Heterotic Eg x E'^ string theory 
(and its dual descriptions) [3, [8]. In such string compactification, moduli fields are not required 
to be at special points, and U(l)-symmetry breaking vev is not assumed to be hierarchically 
small in order to make sure that the trilinear couplings (11.11) are sufficiently small. Thus, this 
solution does not share the unsatisfactory aspect of the solutions with discrete symmetries. 

In section 12.21 a class of compactification of the Heterotic string theory is discussed|§ tech- 
nically, it is to assume that a vector bundle has an extension structure, and various low-energy 

^There are also solutions in the literature where such a U(l) symmetry is broken by a hierarchically small 
expectation value. Such solutions, however, share the same unsatisfactory aspect as those with discrete unbroken 
symmetries. We need to understand why our vacuum is very close to a U(l)-symmetry enhanced point in the 
moduli space. 

^Partial results of sections 12.21 and [2.51 have been obtained in [71 [5] • 
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degrees of freedom are identified with cohomologies of appropriate sub-bundles. We will see 
in this framework that holomorphicity controls mixings between massless states with differ- 
ent U(l) charges. Thus, the selection rule based on U(l)-charge counting is applied for terms 
trilinear in massless states, and the absence of R-parity violating terms (11. ip can be guaran- 
teed from the selection rule. 1-loop amplitudes generate a bilinear R-parity violating mass 
term W 3 fiiLiHu with /ij proportional to supersymmetry breaking (SUSY-breaking), and the 
tree-level contribution can be even smaller. Thus, at the renormalizable level, this framework 
predicts an R parity violation only in the bilinear terms, which is known not to be terribly 
bad in phenomenology [H |TU] . We also find that all the dimension-5 operators that violate R 
parity are generated; the selection rule does not have a predictive power at non-renormalizable 
level. With a theoretical framework controlling all aspects of R parity violation, we can discuss 
how key parameters of short-distance description control the coefficients of various R parity 
violating operators. 

Section [3] is devoted to phenomenology that is expected when both bilinear and dimension-5 
R parity violating operators exist. Although small trilinear R-parity violating couplings can be 
generated in the framework of section [21 it turns out that they are so small that they are rarely 
relevant to phenomenology. Because of negligibly small trilinear R-parity violation, most of 
phenomenological constraints discussed so far are easily satisfied, considerably simplifying phe- 
nomenological study. Remaining constraints from low-energy neutrino masses and washout of 
baryon/lepton asymmetry are briefly discussed in sections [3?T] and [3T2| respectively. Constraints 
on R-parity violating decay of the lightest supersymmetry particle (LSP) are reanalyzed in sec- 
tion 13. 3[ where we exploit the latest understanding of impact of new physics on the Big-Bang 
Nucleosynthesis (BBN). Section IXH is devoted to limits on R-parity violating couplings from 
nucleon decay amplitudes. Although trihnear R-parity violating couplings do not induce too 
rapid a proton decay, squark-exchange diagrams combining bilinear and dimension-5 R-parity 
violating operators still induce nucleon decay. We will obtain a big picture of allowed region of 
parameter space of bilinear-dimension-5 R-parity violation, and find that natural expectation 
of these parameters that follows from the framework in section [2] fits well within the allowed 
region. 

If time is limited, it is best to read only summary sections 12.61 and 13.51 skipping all the 
rest. It is also possible to read sections [2] and [3] separately, as the materials in these sections 
do not require perfect understanding of the contents of the other. 
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2 Theoretical Framework 



Section 12.11 provides basic knowledge in Calabi-Yau compactification of the Heterotic Eg x E'g 
string theory that is necessary later in section [2l We include this mini-review just to make this 
paper self-contained. Thus, there should be no problem in skipping section [TT] and proceeding 
directly to Eil 

2.1 Mini-review on Heterotic String Compactification 

Origin of Gauge Fields and Matter Multiplets 

Effective theories on 3+1 dimensions can be obtained by compactifying the Heterotic Eg, x E'^ 
string theory on a six-dimensional manifold X. A/" = 1 supersymmetry is preserved in a low- 
energy effective theory when X is a Calabi-Yau 3-fold (and hence is a complex manifold). Local 
complex coordinates are denoted by 2" (a = 1, 2, 3). 

The gauge group Eg x Eg can be reduced virtually to any subgroups in low-energy effective 
theories by turning on non-trivial gauge-field background on X. We will refer to the gauge field 
backgrouncifl as V. The gauge group of the effective theory is SU(5)gut for unified theories, 
if the background gauge field configuration is contained within SU(5)' x E'^ C Eg x E^, where 
SU(5)' commutes with SU(5)gut within Eg. 

A super Yang-Mills multiplet on 9+1 dimensions consists of a vector field Am and a gauge 
fermion A vector field on 9+1 dimensions, AM{x,y) (M = O,--- ,9), is decomposed into 
Afj^{x,y) (/i = 0, ■ ■ ■ , 3) and Aa{x,y) (a = 1,2,3); Aa{x,y) is a complex conjugate of Aa 
and is not an independent degree of freedom. Here, (/i = 0,1,2,3) denote four Minkowski 
coordinates, and y"^ (m = 1, ■ ■ ■ , 6) (or simply y) six real local coordinates on X, equivalent of 
three complex coordinates (2;", z") {a = 1, 2, 3). 

Those fields on 9+1 dimensions are decomposed into infinite degrees of freedom on 3+1 
dimensions: their Kaluza-Klein decompositions are given by 

A^{x,y) = J2Aj.^ix)ipi{y), (2.1) 

Aaix,y) = gYM^(l)jix)(pj.aiy), (2.2) 
J 

using mode functions (pi{y) and ip j-aiy)'- Ai-^{x) and (p j(x) are vector and complex scalar fields 
on 3+1 dimensions, respectively. We factored out qyui in (12.21) in setting the normalization of 

■^To be more precise, V means a vector bundle in the fundamental representation that a gauge field back- 
ground in an SU(n) subgroup of defines. 
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the mode functions ^fj-^a] Qym is the gauge couphng constant of the SU(5)gut effective theory 
on 3+1 dimensions. A gauge fermion on 9+1 dimensions, 

v&(x,i/) = (vl>^(x,|/),vl>^(x,|/)), (2.3) 

sphts into a = 1,2,3 part and a = 4 part on a Calabi-Yau 3-fold X, where a and a denote 
doublet indices of left- and right-handed spinors of SO (3, 1), and a label four different entries of 
spinor representations 4+4 of local Lorentz symmetry SO (6). The two parts have Kaluza-Klein 
mode decompositions separately: 

{^T\n=A){x.y) = Y.^hA^)xi{y)r>^''{x)xi{y)). (2-4) 
(vi>^,vi/^)(x,i/) = gYMY.{i^jA^)xj{y).rAx)xjAy)) (a = 1,2,3). (2.5) 

Mode functions eaaXjiv) (being multiplied by a sechsbeir0 Caa) are proportional to those of 
the vector fields fj-a{y) in a supersymmetric compactification, and Aa{x,y) and \E'^(x, ?/) are 
grouped into a Kaluza-Klein tower of chiral multiplets ^i{x, 6, 9): 

A^{x, z, z, e)dz^ = Y,{<Pj + O^J + ---) {x)^j,Mdz^ = J2 ^J^J;adz^- (2.6) 

The remaining (/i = 0, ■ ■ ■ , 3) part and '^'^^ part are also grouped into a tower of vector 
multiplets Vi{x, 9, 9): 

V{x, z, z, 9,9) = J2 {G^^OAj,^ + 9^9Xi + ---){x) ipj{y) = ^ Vnpi. (2.7) 
/ I 

As we assume that gauge-field background is no n- vanishing, mode functions y^j-aiy) and 
fi{y) are not the same everywhere in Eg. The vector field and gauge fermion of the Eq 
super Yang-Mills theory are in the adjoint representation of E^. They split into irreducible 
components 

248 (adj., 1) + (1, adj.) + [(5, 10) + (10, 5)] + h.c. (2.8) 

of SU(5)' X SU(5)gut C Eg. Each irreducible component has its own Kaluza-Klein decompo- 
sition fl2.6l 12. 7p : Mode functions in the (i?', i?)-irreducible component are determined by mode 
equations on the gauge-field background in the R' representation, and hence the spectrum and 

^Hermitian metric haa of a Kahler manifold X is given by haa = Sa=i ^a^aa- 
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decomposition of one irreducible component are different from tliose of another. Therefore we 
use such notations as (i?', R)j or Rj, instead of $j for chiral multiplets. 

On the SU(5)' gauge-field background, massless vector multiplets are found only in the 
(1, adj.) component. Massless chiral multiplets in the SU(5)gut--R representation are in one to 
one correspondence with the zero modes (fij-aiu) in the {R',R) irreducible component. Differ- 
ence between the number of massless chiral multiplets — net chirality — in a Hermitian conjugate 
pair of irreducible components, {R' , R)-{R' , R) , is determined only by topology of the back- 
ground gauge-field configuration V. Phenomenological request is to find a topology of {X, V) 
so that the net chirality in the (5, 10)-(5, 10) sector is three. It then follows that the net 
chirality in the (10, 5)-(10, 5) sector also becomes three; low-enrgy effective theories cannot be 
anomalous if they are obtained by compactifying an anomaly free theory in higher dimensions. 

Vector-like Massless Pair 

Topology of {X, V) determines the net chirality, but the number of massless chiral multiplets 
of each irreducible component can vary for a continuous deformation of gauge field background: 
There can be 3 + m massless chiral multiplets in the 5 representation and m in 5; m can change 
while keeping the net chirality (3 + m) — m = 3. For a given topology of (X, V), however, there 
is still a minimum number of m for the 5-5 sector, and it is not necessarily zeroO See [11] for an 
explicit model, where m > 34 for the 5-5 sector, while there can be no extra massless vector- 
like pair in the 10-10 sector. Similar examples can be found in the literature studying spectra 
on orbifold compactifications. So far, top down principles have been unable to determine the 
topology of {X, V), and hence the minimum number of m's. 

There are a few bottom-up constraints on the number of vector-like massless pairs. First of 
all, there should be at least one vector-like pair in the doublet part of the 5-5 sectoral They are 
identified with the two Higgs doublets of the minimal supersymmetric standard model (MSSM). 
If there were too many SU(5)GUT-charged vector-like pairs, however, they would contribute to 
beta functions of the MSSM gauge coupling constants. Perturbative gauge coupling unification, 
one of the most important motivations of low-energy supersymmetry, would not be maintained 
any more. As long as the "massless" pairs have masses of the order of SUSY-breaking scale or 
higher]^ they have not showed up in experiments, and moderate number of them are tolerable 

^The same thing can happen for the multiplets in the SU(5)gut-10 and 10 representations, in principle. 

^ SU(5)gut symmetry can be broken either by a Wilson line [12] or by a line bundle [13]. Spectra and mode 
functions in a given SU(5)gut representation can be different for different irreducible components of the MSSM 
gauge group. Hence m can be 1 for the doublet part, while m ~ for triplets. We maintain SU(5)gut notations 
in many places in this paper, mainly to avoid cluttered equations, at the cost of sacrificing rigorousness and 
unambiguousness . 

^ Since all the argument based on Calabi-Yau compactification preserves JV — 1 supersymmetry, a vector- like 
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in phenomenology§| 

Superpotential of the Heterotic theory is given by [H] 



W = cMl y ^ A tradj. [acIA - ^^^^^^ , (2-9) 

which is vahd for all the Kaluza-Klein modes $j in (12.61) |15j . Here, c is a numerical constant 
of order unityj^ Mq — 2.4 x 10^^ GeV and 1] is a dimensionless holomorphic 3-form of a Calabi- 
Yau 3-fold X. We treat purely as a background, as we will only discuss what happens within 
the visible sector in this article. 

Kaluza-Klein masses of infinite chiral multiplets originate from d — 2i (A), second derivative 
of (12.91) with respect to A. Writing the superpotential (12. 9p fully in terms of D = 4 chiral 
multiplets, we have bilinear (Kaluza-Klein mass) and trilinear terms: 

W 3 M5(10,5),(TO,5), + Mio(5,10),(5,TO), + Mi(adj.,l),(adj.,l)„ (2.11) 

+ (5, 10) (5, 10) (To, 5) + y" (10, 5) (5, 10) (10, 5), (2.12) 

+ y^(10, 5) (adj., 1) (To, 5) + i/"^(5,T0) (adj., 1) (5, 10), (2.13) 

+ y'"^(adj.,l)(adj.,l)(adj.,l). (2.14) 

Here, multiplets (i?',i?)* represent infinitely many massive chiral multiplets {R',R)i in the R 
representation of SU(5)gut; we will use {R', R)q when we refer specifically to massless modes. 



Those without any * or q stand for both. Chiral 
and (adj., l)o correspond to gauge-field moduli 



multiplets (adj., 1) are SU(5)gut singlets, 
^ Mass matrices such as M5 and Mio are 



of infinite rank. Rank of M5 may be reduced at certain subset of gauge field moduli space. 



pair of multiplets that are massless at a supersymmetric limit may have masses when the supersymmetry is 
weakly broken. 

^Constraints such as FCNC and nucleon decay depend on couplings that those extra particles have with the 
chiral multiplets in the MSSM. Although such constraints can be very severe, we do not discuss them in this 
article. 

^Combined with a Kahler potential, 



K 
Ml 



In {S + S^) - In (^j (T + T^f^ - In (^j 17 A , (2.10) 



this superpotential reproduces a part of the gaugino kinetic term (among other things) with the right dependence 
on a' and Qs- 

Separation between the mass terms and the last two lines is not well-defined; this is because continuous 
deformation of a gauge field background along its moduli space corresponds to changing vev's of massless chiral 
multiplets (adj., l)o, and extra mass terms arise from the last two lines. 
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where extra pairs of chiral multiplets in the SU(5)gut-5 + 5 representations (or possibly in the 
10 + 10) are in the low-energy spectrum [TT] . 
Yukawa Interactions 

The last three lines, f l2.12f[2?T^ . are trilinear interactions involving massless and/or mas- 
sive chiral multiplets, with coupling constants given by overlap integration of relevant mode 
functions. Trilinear interactions involving only massless modes are directly relevant to low- 
energy physics. The first term of fl2.12p contains up-type Yukawa couplings W 3 10 10 if (5), 
and the second term down-type/charged- lepton Yukawa couplings W 3 y'^510H{5). Neutrino 
Yukawa couplings W 3 y'^5NH{5) can only be found in the first term of fl2.13l) . and hence 
the chiral multiplets for right-handed neutrinos are identified with (a subset of) gauge-field 
moduhH (adj., l)o [16]. 

As we have introduced no distinction between C H{5) and L C 5's, however, the second 
term in fl2.12p generically contains the trilinear R-parity violating operators (11. ip as well. Thus, 
a generic {X, V) is not acceptable phenomenologically because of the dimension-4 proton decay 
problem. 

R parity 

The most popular solution to this problem is to impose an R parity. In terms of compacti- 
fication of a super Yang-Mills theory on 9-1-1 dimensions, this is to assume a Z2 symmetry in 
{X, V). We assign odd R-parity for three chiral multiplets 10^ = {Qi,Ui,Ei) {i = 1,2,3) in 
bottom-up model building, but it corresponds to assuming that there are three massless states 
in the (5, 10)~ irreducible component in this context. One further need to assume that there 
are none in the other irreducible components (5, 10) + , (5, lO)"*" and (5, 10)~. Similar assump- 
tions have to be made for irreducible components that are in the 5 and 5 representations of 
SU(5)gut- See Table [U for more details. 

Both moduli (adj.,1)^ and (adj.,l)Q may have non-vanishing vev's without breaking 
SU(5)gut symmetry, but non-vanishing vev's of the latter break the Z2 symmetry. Thus, 
none of (adj., l)o 's should develop a vev in order to maintain an unbroken R-parity. 

With such an unbroken Z2 symmetry of {X, V), down-type/charged-lepton Yukawa cou- 
plings 

W 3 /(lO, 5)0 (5, 10)0 (10, 5)+ (2.15) 
do not vanish, yet the Z2-odd operators 

W 3 A(10, 5)0 (5, 10)0 (10, 5)0 (2.16) 



^^Strictly speaking, right-handed neutrinos do not have to be identified with zero modes. See also the 3+2 
model to be discussed in section \T2\ 
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Table 1: List of chiral multiplets in scenarios with an R parity. In the third column, m = for 
the SU(3)c-triplet parts and m = 1 for the SU(2)i-doublet parts in the third row of this table. 
We listed just the minimum number of chiral multiplets required in effective theory on 3+1 di- 
mensions, ignoring a vector-like (almost) massless pair that can exist without phenomenological 
problems. 

are absent because the overlap integrations for the couplings vanish. 
2.2 4+1 Model and 3+2 Model 

Gauge field configuration V in SU(5)' G Eg on a Calabi-Yau 3-fold X should not be generic, 
since generic configuration gives rise to the R-parity violating trilinear operators (11. ip . leading 
to too rapid proton decay. An R-parity preserving configuration {X, V) with a Z2 symmetry 
is an example of non-generic cases. 

Discrete symmetries other than R parity have also been discussed in the literature as so- 
lutions to the proton decay problem. Phenomenological consequences of these theory can be 
different from those of R-parity preserving ones; the LSP may not be stable, for example. De- 
spite the difference in phenomenology, all the solutions based on discrete symmetries — whether 
preserving R parity or not — are quite similar in philosophy. Enhanced discrete (or continuous) 
symmetries are left unbroken only at special points (or subsets) of moduli space {X, V). There- 
fore, solutions based on discrete symmetries are based on a belief that some dynamics that we 
do not know today will eventually select out vacua with enhanced symmetries, and lift all the 
other part of moduli space. Certainly this belief is not without reason. CP symmetry does not 
have to be preserved in QCD, but once we knew that QCD instanton effects generate a potential 
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of axion, we understood that CP-preserving Ocs. = is the minimum of axion potential. So, 
who can say that history does not repeat itself? 

In this article, however, we neither resort to this belief, nor assume an unknown dynamics 
for vacuum selection. We will present an alternative solution to the dimension-4 proton decay 
problem, which holds at generic points in the gauge-field moduli space. This subsection de- 
scribes a class of compactification that solves the problem. Key ideas were already written in 
[U [S], but a few important clarifications are newly added in this article. We then move on in 
the rest of section [2] to discuss what kind of operators are to be expected in low-energy effective 
theories on 3+1 dimensions. 

4+1 Model 

Let us first suppose that a gauge field background is restricted to an SU(4) x U(l)x subgroup 
of C SU(5)' C Eg. Then, the gauge symmetry is (SU(5)gut x U(1)x)/^5 in an effective theory. 
U(l)x/^5 contains matter parity as a subgroup, Z10/Z5 ~ Z2. Although the matter parity is 
not broken at this moment, we will see later that the restriction on the gauge field background 
is relaxed and the matter parity is broken, yet the trilinear matter parity violating operators 
(jl.ll) are absent. 

Each irreducible component of the i^g-adj. representation in (12. 8p is further split up as the 
group of gauge field background (called structure group) is reduced from SU(5)' to SU(4) x 

U(l)x: 

(5,10) ^ (4,10)-l + (l,10)+^ (2.17) 

(a25,5) ^ (A^4,5)-2 + (4,5)+^ (2.18) 

(adj.,1) (adj., 1)° + (4, 1)"^ + (4,1)+^ + (1,1)°. (2.19) 

Massless modes are identified with various chiral multiplets of the MSSM as in Table [2] in a 
compactification with this class of gauge field background^ C H{5) and Hd C H{5) are 
completely vector like, as in the scenario with an R parity. Massive modes — those that acquire 
masses through W 3 tT{AdA — 2iA {A) A) — have subscripts ^, and massless modes have 0, like 
in Table [H Topology of geometry and gauge field configuration on it, {X, V), should be chosen 
so that the right number of massless modes are obtained. 

^^Here is our naming rule of various chiral multiplets. Since anti-chiral multiplets containing right-handed 
quarks and Icptons arc denoted by U\ and E\ we save N'^ for right-handed neutrinos. Chiral multiplets ^E* 
and '^'^ = 10, 10', 5, N , for example) arise from a Hcrmitian conjugate pair of irreducible components (i?', R) 
and (i?', R) in E^. The same rule is also applied to the 3-1-2 model. In the 4-1-1 model, therefore, H{^) = H{5y, 
but this is not the case in the 3-1-2 model. Since gauge-field moduli (adj.,l)° is a vector-like representation of 
SU(4) X SU(5)gut X U(l)x) there is no distinction between $ and 
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Table 2: List of chiral multiplets in the 4+1 modeL Those with a subscript o are zero modes of 
Kaluza-Klein decomposition, and those with * infinitely many massive modes. Chiral multiplets 
without either one of the subscripts o or * that appear in the text represent both of them. The 
second column shows how the chiral multiplets transform under SU(5)gut D SU(3)c x SU(2) l x 
U(l)y as well as underlying broken symmetries SU(4) x U(l)x iii this model. Only the minimum 
number of chiral multiplets are shown in the third column. The slashes in the third column 
represent the doublet-triplet splitting. For more, see the caption of Table [1] and footnote [181 



Superpotential is obtained by rewriting (12.91) in terms of chiral multiplets in the irreducible 



component 



W = |/fi)1010i/(5)+i/^2)1010'.5: (2.20) 

+ i/f3)io:io:^+yf4)io:io':5 (2.21) 

+ yfi)510iJ(5) + yf2)^(5)10'.iy(5) (2.22) 

+ yf,^5tl0tHi5)+yf,)H{5)10':Hi5) (2.23) 

+ yl,)N5Hi5)+yl,)Nl5lH + y'^,)H{5)<^H{5) + y'^,)5<^5l (2.24) 

+ y;^) ivio', lo: + y[^^Ni 10': lo + y;^) 10$ lo: + y[i) 10: $ loi'^ (2.25) 

+ yl^)N(^N' + yl^)^^(^ (2.26) 

+ MiolOao: + Mio'lO'aO': + MhH.H, + MSM + MmNJK + M$<l>,$,.(2.27) 



An equation (j2.20l - [2?26l) sets up notation for trilinear couplings in the 4+1 model. When we just simply 
use ?/" and y'^ in the text, however, they mean some of j/^-^^.j-^-j and y'^-^^y^^y respectively. Similarly, y^ stands for 
some of y('i)_(4), y(i)_(4) and y'^^y^^^)- 
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Dropping all the interactions involving heavy multiplets, the Yukawa couplings of quarks and 
leptons remain in the superpotential above. On the other hand, trilinear R-parity violating 
operators (11 .ip are not found anywhere in (12.20^2126]) . because 



4,5)+^ ® (4, A^S)"^ ® (4,5)+^ (2.28) 



does not contain a singlet of SU(4) x SU(5)gut x U(1)x- 

The Fayet-Iliopoulos parameter of U(l)x symmetry is given by 



327r2 



Ci(L) A J A J - / Ci(L) cm - ic2(rX) 



s 



vol(X) }^'' ' 2 }^ ' y ' 2 



(2.29) 



which depends on Kahler moduli J (volume moduli) in the tree-level contribution (the first 
term) and dilaton expectation value at the 1-loop level (the second term: a piece well-known 
since 1980's). = 2'K\fa' is the string length. See [I3l[7] for details of the convention. As long 
as moduli fields J and 04 are stabilized by potential other than the D-term of U(l)x symmetry, 
there is no reason to believe that the tree and 1-loop contributions cancel one another. Thus, 
is not expected to vanish. If it is negative, then +5|A^^p in the D-term potential may absorb 
to restore supersymmetry. This scenario is called 4+1 model0 The matter parity Z10/Z5 ^ Z2 
is broken by an expectation value of A^^. 

An order-of-magnitude estimate of U(l)x breaking vev was obtained in [7] . We assume 
that there is no significant cancellation between the tree and 1-loop terms in (12.291) . and that a 
Calabi-Yau manifold X is "isotropic" , that is, its volume moduli are characterized by only one 
typical radius R. Then we can roughly estimate the Fayet-Iliopoulos parameter in terms of 
Kaluza-Klein scale Mkk ^ 1/ R- ■Cx imn determines a vev of A^*^. Canonically normalized 
A^'^ typically develops a vev of order 

which is roughly around the Kaluza-Klein scale [7]. 

The 4+1 model corresponds to turning on a rank-5 vector bundle V that is given by an extension of a 
rank-4 vector bundle 1/4 by a line bundle L: 

(2.30) 

L is a sub-bundle of V , and L ®V and f\^Un are sub-bundles of f\^V and K^V , respectively. 
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Such a large vev of generates Majorana mass terms of right-handed neutrinos A^. Once 
massive chiral multiplets <l>*'s are integrated out, an effective interaction 



u\2 



W 3 ^^^^N'NN'N (2.32) 



is generated [T7]. It is also known [18] that world-sheet instanton effects generate 

W 3 e"^27272727 (2.33) 

in some compactifications of the Heterotic string theory with an unbroken Eq gauge group, and 
(I2.33P contains an interaction of the form 

W 3 —WWNN. (2.34) 

Once the interaction of this form is generated, either from (12.321) or (12.331) . then non- vanishing 
vev of N'^ provides Majorana mass terms of right-handed neutrinos [I9], with masses of order 



Mr = ^ —■. — max 



Vii) 



-T 



\ J 



(2.35) 



Majorana right-handed neutrinos generate small masses of left-handed neutrinos through the 
see-saw mechanism. At the same time, a flat direction SjA^'^p — 5\N\^ + ^;^ = is hfted because 
of the Majorana mass terms. Without a fine-tuning of expectation values of gauge-field moduli, 
we can obtain (A^*^) 7^ while (A^) = 0; (A^'^) 7^ is crucial for neutrino masses as we have 
seen above, and ^A^) = is crucial for the absence of trilinear R-parity violation [7]. 

In the presence of non- vanishing {N'^), trilinear terms of (I2.241[2l26l) involving A^'^ gives rise 
to extra mass terms. We have already seen how the mass matrix in the SU(5)GUT-singlet sector 
is deformed; an interaction (I2.32p with A^'^ replaced by their vev's is regarded as a mass term 
obtained after diagonalizing the mass matrix of $ and A^ [17]. Let us now look at how (A^'^) 7^ 
modifies mass matrices of SU(5)GUT-charged multiplets. 

In principle, mass eigenstates in a given representation of SU(5)gut are mixture of states 
with different U(l)x charges, because the U(l)x symmetry is already broken spontaneously. In 
the SU(5) GUT- (10 +10) sector, the mass matrix becomes 



10: 10- ) (■ 



A/10 






Mio' 




(2.36) 



13 



By this 2x3 matrix, we actually mean an (oo + oo) x (3 + oo + oo) matrix for cases of practical 
interest. U(l)x eigenstates and mass eigenstates are related by a basis transformation^^! 



/ lOo ^ 




10, 




V 10'* ) 





( 



|A/l0'P + lj/'(2)(^')P 





* 




(2.37) 



where 10, are massive mass eigenstates, which consist mainly of linear combination of 10, and 
10'^, and IOq are massless degrees of freedom in the presence of non- vanishing vev of iV . Both 
massive and massless eigenstates 10, and IOq are mixtures of states 10 and 10' with different 
U(l)x charges. 

Mass matrix in the SU(5)gut-(5 + 5) sector is given by 



( Ho H, ) 








Mh 






M5 



H^ 

5o 
V 5, / 



(2.38) 



and one finds that the U(l)x eigenstates contain massless (and massive) eigenstates as in 




1 









f Ho 





* 


* 


II 


5, 





* 






^ 5, 



(2.39) 



A/5 



/\Ah\ 



* * \ 








* * 

1 

* * 



5o 

V 5, / 



(2.40) 



where 5, and 5, are massive mass eigenstates and Hq, 80 and Ho massless degrees of freedom. 
All the massive mass-eigenstates are mixture of states with different U(l)x charges, and so is a 
massless eigenstate H{5)o- But, other massless eigenstates, -f^(5)o and 5o, remain pure U(l)x 
eigenstates, H{5)o and 5o, respectively!^ 



^^Here, we are not very careful in defining the phase of mass eigenstates. 

The disparity between the natures of these massless modes stems from existence of well-defined subbundles 
L, L^Ui and A^U^ that we mentioned in footnote[T4l Massless eigenstates that remain pure U(l)x-eigenstates, 

namely, 5o and Hq, are characterized as H'^{X]L®V) C H^{X\h^V) and H'^{X]A^Ul) C H'^{X; A^V). Other 
massless eigenstates such as IOq or j^o do not have such characterization associated with subbundles [7]. 
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It is important to note that there is a strict rule on the mixing of massless eigenstates. 
Massless eigenstates have their own U(l)x charges in a (N^) hmit. When (A^'^) does 
not vanish, they can pick up interactions of states with different U(l)x charges, only when 
holomorphic insertion of (A^*^) can supply the right U(l)x charge [3, |8]. An H-\ike massless 
eigenstate Hq have a non-vanishing 5 component — (4, 1) entry of the mixing matrix (12.401) — 
because (A^'^) H has the same U(l)x charge as 5. On the other hand, 5-like massless eigenstates 
5o do not have H components — vanishing (2, 2) entry of fl2.40p — because (A^*^) 5 does not have 
the same U(l)x charge as H. Mixing of massless eigenstates IOq is also understood this way. 

Mixing matrices are important, because interactions of mass eigenstates are obtained by 
substituting (12.371) and (12.401) into (I2.20f[2727|) . Since massive mass-eigenstates are generic 
mixtures of states with different U(l)x charges, the U(l)x symmetry is virtually powerless in 
controlling their interactions in the superpotential. Trilinear terms in the superpotential that 
involve only massless eigenstates, however, are still controlled by the U(l)x symmetry, because 
the mixing of massless eigenstates is under the rule above: U(l)x charge can be supplied only 
through holomorphic insertion of {N^)- Three point interactions of massless states do not 
exist, if sums of U(l)x charges of these states in the (A^^) — > limit are positive, because such 
interactions are not U(l)x invariant even after allowing holomorphic insertion of positively 
charged {N'^)- Here, we have a selection rule in the superpotential. 

The R-parity violating trilinear interactions BqIOoSo have positive (+5) U(l)x charges, 
and holomorphic insertion of positively charged (A^^) cannot make them neutral under U(l)x- 
Thus, such interactions do not exist, and the 4+1 model is a solution to the dimension-4 proton 
decay problem. The U(l)x"Charge counting allows an R-parity violating interaction of the form 
W 3 (W) HqWoHo. Indeed, by substituting (12371) and (ICTD into (l2:22|) . we find 

W 3 I ^ '^ =Hol%Ho. (2.41) 
^JM^^ + \y^{N^) P 

But this operator vanishes because of anti-symmetric contraction of SU(5)gut indices and of 
the fact that there is only one down-type Higgs doublet in the MSSM. Therefore, trilinear 
R-parity violating operators are absent in the 4+1 model. 

Finally, let us comment on what it takes in the 4+1 model to have two Higgs doublets 
in the low-energy spectrum, and a //-term in the superpotential. In the doublet parl0 of the 

The doublet and triplet parts of ^^(5) and Hi^) have different KK towers, especially different numbers of 
massless modes, as discussed in footnote IH Mass matrix (|2.38p is intended to be for the doublet part. There 
are no massless modes for the triplet part, so the 1st raw and the 1st and 3rd columns are dropped from the 
mass matrix. 
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SU(5)gut-5 + 5 representations, we need three 5o-like chiral multiplets Lj, one HQ-\ike and 
one Ho-like in the low-energy spectrum. This means that there should be at least 1 and 
m' + 1 zero modes in the doublet part of H{5) and H{5) sectors, and m' and 3 zero modes in 
the 5'^ and 5 sectors (m' > 0}, respectively, m' > is allowed when m' pairs of chiral multiplets 
Hq-Eq acquire large masseo through the second term in (12.241) . Whether m' = or not, we 
need a pair of vector-like massless chiral multiplets in the doublet part of H{5)-H{5). 

The existence of this extra massless vector-like pair is not guaranteed by topology. This 
pair, essentially the two Higgs doublets of the MSSM, can be in the low-energy spectrum for 
two possible reasons. The first possibility is that there is at least one extra pair of massless 
multiplets in the doublet part of H{5)-H{5) sector for generic gauge-field configuration in 
a given topological class of {X, V) . See section 12.11 for more about this case. The second 
possibility is that the background gauge-field configuration of our world is somewhat gpeczaj^ 
and an extra pair of multiplets becomes almost massless for the choice of background. 

It is not hard for the second case for such two Higgs doublets to have a /i-term. Since they 
have trilinear couplings W 3 J2i ^ifiH^HQ with gauge field moduli $/;o (here the lower suffix 
I denotes the /-th modulus), yU-term is generated once the vev's of these moduli fields shift by 
of the order of SUSY-breaking scale. This is essentially the next-to-minimal SUSY Standard 
Model (NMSSM)@ In the first case, such a trilinear coupling is absent (by definition). Instead, 
1-loop diagrams generate such terms as 

This one comes from the first diagram of Fig. [1], and there are also similar contributions from 
loops with multiplets in other representations; see Fig. [H Remembering that holomorphic 
mass parameters in the superpotential have non-vanishing 6*^ components, at least by of order 
M(l -|- 9'^ 771^12) when {W)* /Mq ^ 0, /i-term and i?^-term are generated at 1-loop. This is 
essentially the mechanism in [21] , except that the vector-like massive multiplets in the loop are 
identified with Kaluza-Klein towers here. Note, however, that the /i- and 5^-terms generated in 
this way are known to have a fine-tuning problem [22] . Thus, the problem may be an indication 
that the Higgs sector is a little more complicated than we imagine here, and that a bit of model 
building is necessary as in gauge- and anomaly-mediation scenarios. 



Minimal choice m' = is assumed implicitly in the mass matrix (|2.38p and Table [2l 

Possibly for an anthropic reason [20]. Special choice of gauge- field configuration could follow as a conse- 
quence of moduli stabilization, in principle, but no such claim based on an explicit string construction has been 
made so far. 

^''Hence such moduli are free from moduli problem. 
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Figure 1: Super Feynman diagrams that generate /x-term in the 4+1 model. There are three 
different kinds of graphs; a pair of multiplets in the 10 and 10 representations of SU(5)gut are 
running in the loop (a), those in the 10 and 5 representations of the SU(5)gut in the loop (b) 
(c), and finally, singlets and 5's in the loop (d) (e). 
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r(^T\rt^cit^Y\i' a "l" ion c 






lOo, 10* 
10, 


(1,2,a25)-3 
(1,2, A 5)+^ 


3 

u 


Qi,Ui,E, {i = 1,2,3) 


10'* 

in'" 


(3,1,a25)+2 



n 


- 




(3,2,5)-^ 

fa 2 5)+^ 


3 (+6) 


Di,Li {t = 1,2,3) 


* \"/ 


(1,1, 5)+^ 
(1,1, 5)-^ 


1(2) / 0(3) i+b) 



Hu 


i^o(5), i^*(5) 
^^(5) 


(a23,1,5)+4 
(a23,1,5)-^ 


1(2) / 0(3) i+a) 



Hd 


No, N, 

x 


(3,2,l)-5 
(3,2,l)+5 


1 



heavy RH neutrino 




(adj., 1,1)0 
(1, adj., 1)0 







Table 3: List of fields for 3+2 model. Fields with subscripts o or * can be understood as 
before. The second column shows how the chiral multiplets transform under SU(5)gut as well 
as underlying broken symmetries SU(3) x SU(2) x U(l)5'7 in the 3+2 model. The slashes in the 
third column represent the doublet-triplet splitting, a = 6 = in the minimal choice. We will 
find in section 12.51 that some effective operators can be enhanced depending on whether such a 
non-minimal pair (a 7^ or 6 7^ 0) of (eventually massive) chiral multiplets are in the spectrum 
or not. 



3+2 Model 

There are some variations of the 4+1 model [7], and one of them is called 3+2 model. 
Instead of restricting the structure group of gauge field background to SU(4) x U(l)x C SU(5)' 
at the beginning, one can choose SU(3) x SU(2) x U(l)g7 C SU(5)' as the structure group. This 
restriction on the structure group is relaxed later by U(l)g7-breaking vev, just like in the 4+1 
model. Irreducible components in (12. 8p split into 

(5,a25) ^ (1,2,a25)^^ + (3, 1, A25)+^ (2.43) 
(a25,5) ^ (l,A22,5)-'^ + (3,2,5)-l + (A23,l,5)+^ (2.44) 
(adj.,1) ^ (adj.,l,l)0 + (l,adj.,l)0+(3,2,l)"5 + (3,2,l)+5 (2.45) 

under SU(3) x SU(2) x SU(5)gut x U(l)g7- Each irreducible component has its own towers of 
chiral and vector multiplets. Topology of (X, V) should be arranged so that the number of 
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zero modes of each component is appropriate for a low-energy effective theory; Table [3] shows 
the required number of zero modes as well as notation of chiral multiplets originating from each 
sector. Note that the tower containing the down-type Higgs doublet H{5) is not the Hermitian 
conjugate of that containing the up-type Higgs doublet H{5) in the 3-1-2 model, unlike in the 
4-1-1 model. 

Superpotential of the 3-1-2 model is given by rewriting (12. 9p : 





y^.^lOlOH + y^2)1010',5: + y[^,)10', 10', 


(2.46) 


+ 


yf.^ioiioiH: + 2/^5)10: 10': 5 + yf6)io':io':i5^ 


(2.47) 


+ 


yf,)^105 + yf.^HlO'^H: + y?3)510',5 


(2.48) 


+ 


yf,)H: lo: 8: + yf.^s: w: h + yf.^Ei lo': s: 


(2.49) 


-f- 


2/fi)5iV// + y^.^Hm^ + y'l^^'^'lTHl + y'^.^^W^ 


(2.50) 


+ 


yl,)H<l>H: + yl,-,H^Hl + yl,^f>^K 


(2.51) 


+ 


y;5')¥io:io', +y;^)iv'ioio': +y;^)io$io: +y;^)io',$io': 


(2.52) 


+ 


yl'^^N^N' + y[^^^^^ 


(2.53) 


+ 


MiolO:iO* + Mio'lO'^lO'* + MhHIH. + AfgS^B* + MhKH, 


(2.54) 


+ 




(2.55) 



The trilinear R-parity violating operators (11. ip are absent because 

(3, 2, 5)-^ ® (1, 2, A^S)"^ ® (3, 2, 5)"^ (2.56) 

does not contain a singlet of SU(3) x SU(2) x SU(5)gut x U(l)g7- 

If the Fayet-Iliopoulos parameter of U(l)5'7 is positive, only zero modes Nq, which carry 
negative U(l)5'7 charge, develop non-vanishing expectation values. Majorana mass terms of A^*^ 
are generated in the 3-1-2 model, just like those of are in the 4-1-1 model; see the discussion 
following fl2.32[ [233]) . Small masses of left-handed neutrinos are generated through the double 
see-saw mechanism, since vector-like mass terms of N^^N^ are available from the last term of 
fl2.55p [7]. Super-diagram description of the double see-saw mechanism is the one on the right- 
hand side of Fig. [2l The effective mass scale Mn^ea. that characterizes the low-energy neutrino 
mass through mj, ~ (y'^v)'^ /Mj^^cS. is given by 

M«,eff. ~ ^ (2.57) 

with Mji given in (12.350 . 
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Ho 



Ho 



N N 

; — < — 7^ — > — ( ^(1) 



K Nl 




2/(1) ) ^ X ► A < X ► ( Vd) 



fa) 



A 

(iV) {N) 
(b) 



Figure 2: Diagrams which generate W 3 if (5)iJ(5)55 for the 4+1 model (a) and for the 3+2 
model (b). 



In the presence of non-vanishing vev's of G (3, 2, 1) ^, (oo + oo) x (oo + 3 + oo) mass 
matrix of the SU(5)gut-10 + 10 sector is deformed into 



( 10': lo: ) 



Afio' 




y'ii) W 


A/io 



lo; 

lOo 
10* 



which means that the U{l)q^ eigenstates are related to mass eigenstates as 



10', 
lOo 
10, 

Mass matrix in the 5 + 5 sector becomes 

/ 



* * 

1 

* * 



lOo 
10, 
10, 



( H:^ 5^ 5^ Ho H, ) 



Mr 




> 


yh 


(N) 
(N) 


yh (^) 

vh (N) 








yh 
yh 


W 1/(1) w 
(^> 1/(1) w 





(2.58) 



(2.59) 



(2.60) 



20 



and the U(l)g7 eigenstates have massless eigenstates as components specified by 



Hit 
5o 



5 



(2) 



\ 







Mi 



O 



\ 



M^Mh 



y|M„P + |y^^)(7V)| 



5o 

V 5. y 



(2.61) 



'\MH?+\y'(^){N)\ 

Using (I2.59P and (I2.6ip in (12.480 . one finds that the trilinear R-parity violating operators (II. ip 
involving massless eigenstates are not generated. The essence is that U(l)g7-eigenstates H{f>) 
do not have the massless eigenstates 5o as a component, and U(l)g7-eigenstates 10' do not 
contain massless eigenstates IOqO Eigenvectors for massless eigenstates in (I2.59[ I2.6ip have 
non- vanishing entries only when holomorphic insertion of U(l)g7-breaking vev of can fill the 
gap of U(l)g7-charges, just like in the 4+1 model. Therefore, trilinear interactions involving 
only massless eigenstates are subject to the selection rule based on holomorphic insertion of 
(A^) and U(l)5'7-charge counting. 



In the doublet partlfj of SU(5)gut-5 + 5, one Hu C Hq, one Hd C Hq, and three Li C 5o 
should remain in the low-energy spectrum, although and {Hd-, Li} are in a pair of vector-like 
representations of the MSSM gauge group. Thus, the mass matrix (I2.60p somehow has to have 
a reduced rank for the doublet part. Note that it is actually an ((oo + a + oo + (l + 6) + 
oo) X ((1 + a) + oo + (3 + 6) + oo + oo)) matrix for the doublet part@ although it looks like 
5x5; here a and h denote the number of massless modes in ^'^-H and if-5 vector-like pairs, 
respectively, that become massive after spontaneous ^7(1)5-7 breaking. It is supposed to have 
rank oo + a + 00 + 6 + 00, not 00 + a + 00 + (1 + 6) + 00. 

First of all, in order for three Lj's to remain in the low-energy spectrum, the (1 + 6) x (3 + 6) 
submatrix — (4th, 3rd) block — of (I2.60p should have rank 6, not (1 + 6). Such rank reduction 



"^^ The 3+2 model corresponds to a rank-5 vector bundle V given by an extension of a rank-3 bundle by 
a rank-2 bundle C/2: 

^ C/2 ^ V" ^ C/3 0. (2.62) 

U2 is a sub-bundle of V , and U2 ®V that of l\^V . The massless eigenstates IOq and 5o are characterized as 
H'^{X\ C/2) C H^{X] V) and H^{X] U2®V) C H^{X; ^'^V), respectively, and have restricted interactions than 
general massless states in H^{X; V) and H^{X] A^F), respectively. 
See footnotes ISl and [T7I for discussions of doublet-triplet splitting. 
^•^For the triplet part, (1 + 6) and (1 + a) are replaced by h and a, respectively. Rank reduction is not necessary. 
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of mass matrix can happen for a generic point on the gauge-field moduh space, as we have 
reviewed in section 12.11 and remarked as the first possibihty for the hght two Higgs doublets 
in the 4+1 model. This means that there are one linear combination of Hq^s and three of Sq's 
that do not have couplings of the form W 3 y'^ HoNo5o. 

Secondly, another massless doublet Hii C Hq also has to remain in the low-energy spectrum. 
In general, i^o-type zero modes and -f^o-type zero modes effectively have mass terms of the 
form W 3 [(2/(1) (^o))(2/(2) {^o))/M^]HoHo, and become massive. It may be, however, that 
this mass matrix has a reduced rank for generic moduli field (including A^'s) value, just like we 
assumed above. Alternatively, the effective mass parameter happens to be small at the value 
of (iV^'s where moduli are stabilized. While the first possibility is aesthetically better, the 
origin of fi parameter is not explained because the zero modes Hq and Hq are not coupled to 
moduli fields that might play a role of the singlet field of the NMSSM. The second scenario 
seems to involve a fine-tuning to get the /i parameter small0 some anthropic selection may be 
responsible for it [20] . 

Remark 

There are some other variations of the 4+1 and 3+2 models [7], but they are essentially the 
same in the mechanism of eliminating the trilinear R-parity violating operators. Thus, we use 
these two models as illustrative examples representing a theoretical framework that can be an 
alternative to R parity. 

2.3 Dimension-5 Proton Decay Operator 

The idea in section [2l2l was that the mixing of massless modes is governed by holomorphic inser- 
tion of vacuum expectation value of chiral multiplets. The trilinear R-parity violating operators 
are absent and the dimension-4 proton decay problem is solved, without an unbroken discrete 
symmetry such as R-parity. Supersymmetric extensions of the standard model, however, has 
another problem. If the effective superpotential contains 

AW= -^QQQL + ^EUUD = 0i + 02 C -^1010105, (2.63) 

these terms violate both baryon number and lepton number, and hence proton decays. Thus, 
we discuss in this subsection whether or not the theoretical framework in section 12.21 predicts 
the operators above in the low-energy effective superpotential. 

^''There may be multiple iV's that acquire vev's, and there are infinitely many Kaluza-Klein states 
Thus, the effective /i-parameter being small does not necessarily mean that y'^ (^N'j itself is small. 
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Figure 3: Dimension-5 proton decay operators in W 3 1010105 are generated by a diagram 
(a) in the 4+1 model. Although they appear to be generated in the 3+2 model as well through 
the diagram in (b), it turns out that they are actually not. See the text for explanation. 



4+1 model 

In the minimal supersymmetric SU(5) GUT model, the dimension-5 operators above are 
generated by integrating out colored Higgs multiplets. The essence is that the pair of colored 
Higgs multiplets in the SU(3)c-3 + 3 representations has a vector-like mass term. The 4+1 
model in section 12.21 also shares the same property. The up-type and down-type Higgs doublets 
of the MSSM originate from irreducible components H and H in the SU(5)gut-5 + 5 repre- 
sentations, and they form a Hermitian conjugate pair even in the Eg gauge group. Therefore, 
the superdiagram in Figure El (a) generates the term AW = 1010105 in low-energy effec- 
tive theory, similarly to the minimal supersymmetric SU(5) GUT model. In particular, the 
QQQL + EUUD part relevant to proton decay is generated by colored Higgsino exchange 

There are a couple of differences, however, between the minimal supersymmetric SU(5) 
GUT model and the 4+1 model. In the 4+1 model, propagating between the two vertices 
of Figure [3] (a) are not just one pair of colored Higgs multiplets. All of massive Kaluza- 
Klein excitations contribute to the amplitude. These infinite number of contributions add up 
coherently. Thus, we have an effective operator 

Weff. 3 V ^^lOolOolOoSo = -^lOolOolOoBo, (2.64) 

with massless eigenstates in the external lines. As this effective operator is generated in the 
4+1 model without a mixing induced by {n'^), we do not make an explicit distinction between 
the U(l)x-eigenstates and mass-eigenstates. Yukawa couplings yi are determined by overlap 
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Scalar colored Higgs exchange is also relevant to QQQL + EUUD proton decay. 
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integration of mode functions of three relevant states; 

= ^^YY / '^'y^' (X07™^0;r.X/) , (2.65) 
VOl(_A j Jx 

where xo and (po-m are wavef unctions of massless modes to be used in the external lines, and vr 



that of /-th massive Higgs multiplet in the internal line with mass Mh-i- The overall factoi 
Qyu originates from the normalization convention of the mode functions in (12.21) . 

Sum of infinite contributions in fl2.64p is treated better in an equivalent description 
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Meff VOI(X) 



^YM 

vol(X) 



Jx vol(X) Jx ^ Mh-i 

^ ' d'y\^o;nrXo){y). (2.67) 



dMX07"'^0;m){y) 



X 



(27r)6 ^ 



X 



Discrete summation labeled by I in the first line is approximated by continuous momentum 
integration in the second line. This is a reasonable thing to do, as long as we are interested 
in contributions from highly excited Kaluza-Klein states. The factor in the square bracket in 
fl2.67p is nothing but Green function over the internal space X, and its short- distance singularity 
7^(9fc {l/\y — y'\^) corresponds to summing up contributions from infinitely many Kaluza-Klein 
states. For zero-mode wavefunctions Xo and fo;m without particularly singular or rapidly 
varying behavior, integration over \y — y'\ in fl2.67l) is dominated by the long-distance region, 
\y — y'\ ^ i?, not by the short-distance region. Thus, l/Mefj. is finite. We also learn here 
that contributions from low-lying triplet Higgsino exchange dominate the amplitude 1 / Mcs_ of 
dimension-5 proton decay operators. 

Experimental limits on the dimension-5 operators fl2.63p are set by several decay modes of 
proton. t{p + v) ^ 2.3 x lO'^^yrs. [23], one of the most important ones in the minimal 

supersymmetric SU(5)gut niodel, roughly corresponds to 



Meff. > 10'^ GeV. (2.68) 

Limits from other decay modes are somewhat different, but not by several orders of magnitude. 
All kinds of GUT models are marginally in conflict with this constraint [2U [25] as long as the 
two Higgs triplets are vector like in any underlying symmetry groups. This property is shared 
also by the 4-1-1 model If l/Meg, is approximated by contributions from a few lightest triplet 



Mode functions should be normalized so that the kinetic terms of chiral multiplets are normalized canoni- 
cally. The expression (|2.65p contains e<^>/*^G, where {K) /Mq is a vev of (|2.10p . 
^''The R-parity preserving scenario explained at the end of section [2?T1 also does. 
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massive Higgsinos in the Kaluza-Klein tower, l/M^s, ~ y'^yj / My^]^, the experimental hmit is 
translated to 

This constraint seems a little severe if the Kaluza-Klein scale Mkk is around the GUT scale 
Mqut ~ lO^^GeV, and the trilinear couplings are of the order of the Yukawa couplings of the 
second generation quarks, yc ~ 10"^ and y^ ~ 10"'^. However, y^'^^'s in f l2.64[ I2.69P in the 4+1 
model are not related to the Yukawa couplings of quarks and leptons by SU(5)gut symmetry, 
y/'s are calculated in (12.651) by using mode functions of Kaluza-Klein triplet Higgsinos Xiiv)) 
whereas the observable Yukawa couplings are calculated by using mode functions of two massless 
doublet Higgsinos. Since y/'s are "Fourier transform" of zero-mode wavefunctions (xo7"*V^O;m) 
on a compact internal manifold X, they become very small for higher Kaluza-Klein states. 
Therefore, it will not be conservative to exclude the 4+1 model with Mkk ~ Mqut- 
Remark 

There is a theoretical side remark here, before moving over to the 3+2 model. In the process 
going from (12. 9p and/or (I2.20f[2l27|) to (12.641) . we are rewriting IPI effective action: from the one 
written in terms of (Kaluza-Klein decomposition of) Yang-Mills multiplets on 9+1 dimensions 
to the one in terms only of Kaluza-Klein zero modes on 3+1 dimensions. The dimension-5 
operator (I2.64p is necessary as a IPI vertex when IPI effective action is written in terms only 
of zero modes, although it was absent in (12.91) where all Kaluza-Klein modes appear in IPI 
effective action. It is therefore almost trivial that the sum of infinite contributions in (I2.64p is 
finite; the amplitude described by (12.641) and by Figure [3] (a) is nothing more than a tree-level 
scattering amplitude of super Yang-Mills theory with the external states having wavepackets 
of Kaluza-Klein zero modes. 

3+2 model 

Let us now study whether the dimension-5 proton decay operators (12.631) are generated in 
the 3+2 model. In this model, U(l)q7-eigenstates H and H are not Hermitian conjugate, and 
hence the dimension-5 proton decay operators are absent in the effective theory in the (A^) = 
limit. Thus, whether those operators exist in low-energy effective theory depends on how the 
vev's of can be inserted in the effective superpotential. 

U(l)q7-breaking vev's <^A^) induce mixing between states with different U(l)g7 charges. The 
mixing can take place among massive states propagating in the internal line as in the diagram 
in Figure [3] (b). On the other hand, the massless eigenstates IOq in the external lines remain 
pure U(l)g7-eigenstates, as we saw in (I2.59p . Although iJ^-type U(l)g7-eigenstates contain the 
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massless eigenstates 5o as components, there is no interaction involving 10 and simultane- 
ously in fl2.48p . Thus, mixing in the internal line is responsible, if the dimension-5 proton decay 
operators are to be generated. 

Let M/j denote the mass matrix in (12.601) . and U and V unitary matrices that diagonalize 
Mij- 

Urn Mij Vjk = ShrMk • (2. 70) 

Then the dimension-5 proton decay operators in the effective theory are given by 

^^^ 3 J2 Y^ i^y'jyjKmUi,!) ^^^^^^^^^-^^ ^ -^IOoIOoIOoBo; (2.71) 

/GH(5),Je/?(5) K ^Wcff. 

the prime in the summation means that the label K runs only over massive states. The effective 
coupling 1/Meff. can be expressed in a much simpler way. 

Vi) yHi,)KM-K'UKHi.)l yh = yfi) (^"')/?(5)H(5) Vh- (2-72) 



Thus, although the unitary mixing matrices U and V are not given by rational functions of 
parameters in the superpotential, the coupling l/Mefj. is in the effective superpotential. 

Let us examine the matrix M~^, first with the simplest case a = 6 = in Table [31 This 
corresponds to a case where there is no pair of chiral multiplets in the SU(5)gut-5 + 5 rep- 
resentations that acquire masses after (iV)'s become non-zero. First, we are concerned only 
about the 3 + 3 part of the mass matrix in the 5 + 5 sector, because the dimension-5 proton 
decay operators are generated only through triplet exchange diagrams. Second, three massless 
chiral multiplets Di C 5o are taken out of the triplet part of the mass matrix (12.601] . because 
only massive states are integrated out. Now, the mass matrix (12.60p has become 

/ Ms_ \ 

Mu = y^2) (N) M5_ , (2.73) 
V y'(,^{N) Mh I 

and we treat it as if it were a 3 x 3 matrix. The index I runs over {ff'^(3), 3^, if(3)} C 
{i^'=(5),5^i7(5)}, and J over {i7(3), 3, i/^(3)} C {i/(5), 5, i/^(5)}. There is no qualitative 
problem in approximating the summation over contributions from infinite Kaluza-Klein states 
by a sum over three distinct Kaluza-Klein towers, as we have discussed before for the 4+1 
model. The entire contribution from a tower is finite, and only those from low-lying massive 
states have sizable contributions. It is now easy to see that the cofactor of (-ff(3), -^(3)) element 
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of the matrix fl2.73p vanishes, and so is {M~^)ff(^^^ff(^^y Thus, the dimension-5 proton decay 
operators are not generated in the 3+2 model with spectra characterized hj a = b = 0. 

Individual contribution from a given mass-eigenstate in fl2.7ip . VjkM^^Uki without sum- 
mation over K, does not vanish, as the Feynman diagram in Figure [3] (b) indicates. Some of 
those contributions are of order M{y'^ (A^))*^/|M|'^, where M Mfj^^ jj. But, all these contri- 
butions cancel, and the total amplitude l/Mefj. vanishes. The total l/Mcfj. should be given by 
a rational function of parameters in the original superpotential, as we saw in fl2.72p . Vev's of 
anti-chiral multiplets {y'^ (^N^y^ should not survive cancellation. 

Similar study can be carried out for cases with either/both a or/and b is/are non-zero@ 
The 5x5 matrix fl2.60p without rank reduction is used instead of (12.730 for the case with 
a 7^ and 6 7^ 0, and 4x4 matrices are used for the two other cases. The cofactor of the 
(if(3), ff(3)) element turn out to vanish for all these matrices, and hence the dimension-5 
proton decay operators are not generated for any of these cases. Thus, the dimension-5 proton 
decay operators are absent in low-energy effective superpotential of the 3+2 model, independent 
of whether 5'^-ff(5) like and/or i/(5)-5 like states with masses of order y'^ (^N^ exist at high 
energy or not. 

The absence of the operator in the 3+2 model can be understood in terms of U(l)g7-charge 
counting. Requiring that the effective coupling 1/Meff. is a rational function of parameters and 
vev's of the original superpotential, and that the effective superpotential is invariant under 
the U(l)g7 symmetry, the only possible form is l/Meg. ~ y^y'^M / [y'^ (A^))^, where M is some 
mass scale. Unless there are massive states in the SU(5)gut-5 + 5 representations with their 
masses given by ~ O^i^y" (iV))^/M), such terms are not generated in the low-energy effective 
superpotential. This could have been the case if there were extra pairs of zero modes of H{6)- 
type and -ff(5)-type, in addition to those listed in Table [3l We have ignored such a possibility 
so far, because the dimension-5 proton decay operators with 1/Mcs. ~ U^U'^M^/ {y^ (^))^ seem 
too large compared with the experimental limits, if Mg ~ Mkk ~ ^gut and y^ (^N^ <^ M^k, 
despite large uncertainties associated with trilinear couplings y^ and y'^ that involve Kaluza- 
Klein states. 

2.4 Bilinear R-parity Violation 

Let us now study consequences of R-parity violation in the framework presented in section 
12. 2[ We begin with a study of mixing between Li C So and Hd C H{5)o', they are in the 

^^In these cases, extra massless multiplets in the SU(5)gut-5 + 5 representations appear in the (iV) 
limit. 
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(a) 



10 

(b) 



Figure 4: Diagrams which produce K 3 5^^ H for the 4+1 model (a) and for the 3+2 model (b). 

same representation of the MSSM gauge group S\J{2)l x U(l)y C SU(5)gut, and have distinct 
symmetry charges only in the presence of an unbroken R-parity or U(l) symmetry. Any inter- 
actions involving R-parity violation (and hence a vev of A^'^ (in the 4+1 model) or (in the 
3+2 model)) have a chance to induce mixing between Lj's and H^. 

We have already seen in section 12.21 that the vev's of A^'^ [resp. A^] deform mass matrices 
in the superpotential in the 4+1 [resp. 3+2] model. Consequences of the deformed mass 
matrices, however, were quite limited in the superpotential at the renormalizable level. Mixing 
of massless eigenstates are under strict control of the spontaneously broken U(l)x [resp. 11(1)5-7] 
symmetry, and the trilinear R-parity violating operators (11.11) are not generated for massless 
modes. Although the U(l) symmetry does not rule out the other trilinear R-parity violating 
operator of the form fl2.4ip . it is absent in the MSSM because there is only one down- type 
Higgs doublet. Yukawa couplings of quarks and leptons will have different values from those in 
fl2.20f[2l24|) [resp. fl2.46f[2750|) ]. because the massless eigenstates may not be exactly the same 
as the original U(l)-eigen zero modes. Those are all the consequences. 

Vev's of A^'^ or A^, however, generate mixing in the Kdhler potential as well. Kinetic mixing 

fsTeff. 3 e5^if(5) + h.c. (2.74) 

is generated in both the 4+1 and 3+2 model at 1-loop by super Feynman diagrams in Figure HI 
Those diagrams give rise to operators. 



167r2 |M|^ 



4+1 model i^eff.-^ / 5t^ + h.c., (2.75) 



3+2 model J^es. ~ ig^2 ]m\^ + ^.c. (2.76) 
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|Mps in the denominators stand for the largest among jMgp, |Miop and |Mjjp, because that 
is where the dominant contribution comes from in loop momentum integration. There are two 
other kinds of contributions as well, but they are quite similar to those above; see the caption 
of Figure [H We treat 1-loop amplitudes here, as if only finite number of massive particles ran 
in the loop. 

It is true in the context of string compactification that infinitely many Kaluza-Klein particles 
and stringy states are also running in the loop. The argument above estimates contributions 
only from low-lying Kaluza-Klein multiplets. Unless the remaining UV contributions exactly 
cancels the IR contributions above, however, the total 1-loop amplitude does not vanish. Since 
the IR and UV contributions are likely to depend on geometry differently, it is unlikely that 
they cancel. Thus, as long as there are non-vanishing IR contributions, it is likely that the 
total amplitude does not vanish. This is what we can guess from the argument using Feynman 
diagrams. 

Although the UV contributions to e have not been discussed yet, it is sufficient to treat it 
symbolically in seeing that its effects disappear from renormalizable interactions of the massless 
modes of the MSSM. By a non-unitary basis transformation 

(f)^(-^O(f)' 

the kinetic terms are diagonalized, yet the R-parity violating trilinear operators f 1 1.1 1) are absent 
when the superpotential is rewritten in terms of newly defined chiral multiplets H' and 5'. 
Although operators of the form (12.411) now exist when written in the new chiral multiplets, 
there is no such term that consists only of massless multiplets in the MSSM because of anti- 
symmetric contraction of SU(5)gut indices [8]. Jacobian of this field redefinition in path integral 
is trivial, and there are no other effects from this redefinition. 

Similarly, 1-loop amplitudes generate kinetic mixing among chiral multiplets in the SU(5)gut- 
lO-flO representations. An appropriate field redefinition, however, can diagonalize the kinetic 
terms without generating R-parity violating trilinear operators that involve only massless mul- 
tiplets in the MSSM. In the 3+2 model, for example, an effective Kahler potential term 

~ e 10^10' + h.c. (2.78) 

is generated, with a coefficient e proportional to (iV). A new basis (10', 10)' is chosen by 




(2.79) 
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(a) (b) 

Figure 5: (a) is a typical diagram for (12.811) in 4+1 model, (b) is a typical diagram for (I2.82p 
in 3+2 model. 

so that the kinetic terms become diagonal. Note that the redefinitions of chiral multiplets (I2.77P 
and (I2.79P use parameters e, which involve only holomorphic vev's of either or A^. This is 
why the basis transformation matrices (I2.77P and (12.790 are lower triangular. It thus follows 
that the mass matrices (12.381 12.60p and (I2.58P remain lower triangular, and whole argument we 
have had so far in sections 12.21 and 12.31 does not have to be changed qualitatively§ 
The Kahler potential of the effective theory may also have a bilinear term 

K^s,3 c5H (5) +h.c. (2.80) 

1-loop super Feynman diagrams in Figure [5] generate 

4+1 model i^eff. 3 ' \J. '° 5iJ + h.c., 2.81 

3+2 model K,s. 3 + h.c, (2.82) 

where |M|'s in the denominator are the scale where the dominant contribution of the loop 
momentum integration comes from field theory on 3+1 dimensions. Similar contributions come 
from two other kinds of diagrams with particles in different SU(5)gut representations running 
in the loop, just like several kinds of diagrams in Figure [TJ In both models, the coefficient c in 
(I2.80p is roughly of the form 

lovr^ M 



^^We will make an order of magnitude estimate of e at the end of this section, and argue that the additional 
deformation to mass matrices is quantitatively unimportant compared with the tree-level deformation discussed 
in section 12.21 
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where phases are ignored, and (A^) means {N^) in the 4+1 model and ^A^) in the 3+2 modeL 
The estimate (]2.83p takes account only of 1-loop amplitudes with low-lying Kaluza-Klein mul- 
tiplets in the loop, and there will also be UV contributions where higher Kaluza-Klein states 
and stringy states are in the loop. Although we have not estimated the UV contributions, it is 
unlikely that the UV and IR contributions cancel almost exactly, as we discussed for the kinetic 
mixing K 3 5^ H. Thus, it is likely that c does not vanish in low-energy effective theory. 

Holomorphic mass parameters M are expected to have 6"^ component in a vacuum with 
broken supersymmetry, at least of order M(l + 9'^m^j2)- Therefore, the anti-holomorphic mass 
parameters in fl2.81l 12.821) have 9"^ components, and hence the 1-loop amplitudes generate an 
R-parity violating bilinear term 

W^eff. 9 cm3/25i/(5). (2.84) 

In the SU(3)c triplet part, this additional mass term AVT = /ijZ)jiir(3)* is not signifi- 
cant. It deforms the mass matrix including W 3 MhH'^{^)*H{2>)^, in the SU(3)c-3 + 3 sector, 
and induces mixing between H^{Z)^ and D. Massless eigenstates Di pick up interactions of 
H'^{2)^, C H'^{5)^: with mixing coefficients —fii/Mn- Thus, the massless eigenstates have trilin- 
ear interactions 

^ -vUitn^'UD - vU-^^b,QL, (2.86) 

in the 4+1 model, where H^{f>) = H{^)- When ~ Mkk is around the GUT scale, effective 
R-parity violating coupling A" ~ —y'^Hi/MH is very small, and irrelevant to phenomenology 
except in nucleon decay processes discussed in section 13. 4[ This coupling is not even generated 
in the 3+2 model, because H'^ is different from H. 

In the doublet part, however, the additional mass terms, 

Wes. 3 Cim3/2LiHu = liiLiHu (2.86) 

with Cj roughly given by (12.830 . are important in phenomenology. Since Hu only has a small 
mass term AIV = ^qH^Hu-, with the /i-parameter /iq of the order of the electroweak scale, the 
mixing angle of Li-Hd mixing, ~ /Uj/yUo, is not so small as fii/Mn- Although this mixing angle 
is quite smaller than 0{1) because of the suppression factor c in (12.831) . it is much larger than 
fii/Mn and can be phenomenologically significant. 

In addition to Li-H^ mixing in massless states, Li have mixing of the order fii/Mn with 
massive states in the doublet part of H'^{5)^. In the 4+1 model, this mixing generates trilinear 
R-parity violating interactions 

P^eff. 3 -yt,)-^UEL - yf^^J^UQD, (2.87) 
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since H^{5) = H{5) in the 4+1 model. If Mh ~ A^kk is around the GUT scale, however, 
effective trilinear R-parity violating couplings A ~ A' ~ —y'^^^'^^i/MH are so small that it is 
negligible in phenomenology compared with fii/fiQ contribution from fl2.86p . Thus, even in the 
4+1 model, trilinear R-parity violating couplings fl2.85[ [2^87|) are too small to have phenomeno- 
logical significance, and the bilinear mass term (12.861) is virtually the only phenomeno logically 
relevant R-parity violation at renormalizable level |fj 

We have seen so far that the trilinear R-parity violating operators fll.ip are either absent 
[in the 3+2 model] or negligibly small[in the 4+1 model] in this framework. This framework 
predicts bilinear-dominated R-parity violation at renormalizable level, which was assumed in 
[To] . The dimension-4 proton decay problem is absent in bilinear-dominated R-parity violation, 
as baryon number symmetry is preserved (apart from a small breaking in fl2.85p ). So, it is easier 
to satisfy other phenomenological constraints as well in bilinear-term domination, comparing to 
the case of trilinear-term domination; as long as yUi/yu's (and bilinear R-parity violating terms 
in the SUSY-breaking potential) remain sufficiently small, everything is fine. The only question 
is why /ij//i's are small. 

In our framework, coefficients Cj's are all suppressed by {y'^ {N))/M as in [SHI], /ij = Cim^i2 
becomes even smaller for smaller gravitino mass |6] Even if gravitino mass is not very small, 
Cj's involve an extra 1-loop factor in (12.830 since those operators are generated at 1-loop. Cj's 
may be suppressed further by some ratio of mass parameters that we find in (12.811 12.820 . 
Therefore, there are many reasons for /ij to be small in this framework. 

As we have made it clear, however, the rough estimate of the 1-loop amplitude (I2.83P only 
accounts for a partial contribution. To see if we can rely on the estimate, let us first discuss 
if there are not any tree-level contributions at all. Secondly, we also need to see if the infinite 
Kaluza-Klein particles and stringy states in the 1-loop amplitude could give significantly larger 
contributions than (12.830 . 

In order to find out whether there is a tree-level contribution, it is desirable, in principle, to 
calculate a sphere amplitude on a Calabi-Yau background. In reality, though, such a calculation 
is rarely available except in orbifold limits of Calabi-Yau 3-folds. However, tree-level IPI 
effective action has been calculated for the Heterotic x E'^ string theory on a fiat 9+1 
dimensional spacetime in a'-expansion. Here, a' = and Mg is the string scale. If fiiLiHu 

mass terms are obtained from dimensional reduction of the tree- level IPI effective action on 

■^°In section [5?^ however, we will see that the operator (|2.85p can be important for some choice of parameters. 

■^^Possible contributions from messenger sector fields in gauge mediation models have to be studied separately. 
This issue is not covered in this article, since such contributions will depend very much on details of models of 
gauge mediation. 
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9+1 dimensions, it is likely that they are obtained also in sphere amplitudes for compactified 
models. If /ij's obtained from dimensional reduction are proportional to a positive power of M^, 
it is likely that they are, too, in compactified models, because we expect that the process of 
compactification introduces only Mkk- dependence, not extra Mg dependence. Thus, we content 
ourselves with guessing whether yUj's are generated at tree- level, and if they are generated, how 
they depend on Ms, by using dimensional reduction of IPI effective action on the flat 9+1 
dimensional spacetime. 

In order to find out whether dimensional reduction gives rise to /ij's proportional to gravitino 
mass, we need to assume an origin of non- vanishing (W*). It will depend on the assumption 
whether there are tree-level contributions or not. In this article, we study only one possibil- 
ity for the origin of (W*), just to illustrate what one should consider to guess the tree-level 
contributions to /ij's. 

Let us suppose that (W*) originates from a 3-form flux J^Q AH [26]. If a mass term of 
left-handed chiral fermions is to come through dimensional reduction, then there should be two 
gauge fermions = Yla=i 2 3 ^"^sa ^ term of effective action on 9+1 dimensions. Since the 
mass parameter is supposed to be proportional to gravitino mass, a 3-form vev {Hap^) should 

cf 

also be involved. In the 4+1 model, we further need a vev of an anti-chiral multiplet , 
which originates from a vector field vev. In order to write down such a term that is at least 
invariant under the holonomy SU(3), it must contain at least three derivatives; candidates of 
such terms that appear first in a'-expansion are of the form 

^ ii {D^^<s,Dp^<pD^A^)Hse^d^^y, (2.88) 



with SU(3) indices contracted by metric haai holomorphic 3-form Vtap-y and its Hermitian 
conjugate ^ap^- This term is proportional to l/gl because it is suppose to be in the tree-level 
effective action. A dimensionful coefficient \/a' was supplied based on dimensional analysis. 

One can see, however, that such terms cannot be made S0(6) ~ SU(4) invariant. To see 
this, note first that the SU(3) indices in fl2.88p can be contracted by one fij^^ and some metrics 
haa- Since SU(4) contains SU(3) x U(l), (12.881) should be neutral under the U(l), if it is to be 
made invariant under SO (6) ~ SU(4). Remembering how spinor and vector objects transform^ 
under SU(3) x U(l), one finds that (12.881) with an VL^^^ has +6 units of the U{1) charge. Thus, 
(I2.88P cannot be made invariant under SU(3) x U(l) C SU(4), no matter how the indices are 



A spinor field ^'a is in the 3^ representation of SU(3) x 11(1), vector-type objects such as and 
transform as 3+^, and anti- vector- type objects such as Da as 3~^. 



33 



contracted by metric. Needless to say, there is no term of the form fl2.88l) that is invariant 
under S0(9,l). 

There is no tree-level contribution even from higher order terms in O^a'D"^) expansion. This 
is because DaD^ is neutral under U(l), and so is {DaY accompanied by fia/j^. 

In the 3+2 model, on the other hand, it is vev's of chiral multiplets that are inserted in 
f l2.82p . Thus, a Lorentz invariant operator can exist in 9+1 dimensions at tree level, 

^tr(^,(r%)"'v^,), (2.89) 

which is nothing but a part of the kinetic term of gauginos in 9+1 dimensions. This term 
has become a part of the superpotential (12. 9p . But it was a part of the assumptions of the 
3+2 model that the massless doublet -f^o remains in the low-energy spectrum because of rank 
reduction in the coupling W 3 u'^HqNq^q. Therefore, the existence of (12.891) in 9+1 dimensions 
is irrelevant to the question of whether the Li-Hu mass term is generated at tree level; even 
if the vev of gauge field moduli A^o's shift by of order 171^/2 in the presence of SUSY-breaking, 
nothing happens because the coupling is absent. Higher order terms in a'-expansion may give 
rise to the tree-level 5o--f^o mass term in the 3+2 model as well; 

tr (D^'^^DfsA^D^^^) Hs^c (2-90) 

can be made invariant under the holonomy SU(3) after the indices are contracted by haa, 
and ^af3^- Following the same argument as in the 4+1 model, however, there is no way making 
(I2.90p invariant under SU(3) x U(l) C SU(4) 2^ S0(6). There are no higher order terms in 
0{a' D"^) expansion in the S0(9,l)-invariant effective action whose dimensional reduction gives 
rise to the tree-level contribution, just like in the 4+1 model. 

The other issue is the contributions to the 1-loop amplitude with infinite Kaluza-Klein 
particles and stringy states in the loop. The 1-loop amplitudes (12.811 12.820 are obtained by 
treating Kaluza-Klein towers as if they contained one (or a finite number of) Kaluza-Klein 
excitation(s) in each tower. Such amplitudes are UV-finite, because a superficial degree of 
divergence is —2 and —4 in the 4+1 and 3+2 model, respectively. But Kaluza-Klein towers 
contain infinite degrees of freedom on 3+1 dimensions, and their contributions are approximated 
by integrating loop momentum in lOD, not in 4D. Thus, the superficial degree of divergence 
becomes +4 and +2, respectively. These UV divergence is tamed and made finite in string 
theory. The most naive guess is that the estimate (12.831) is multiplied by a factor of (M^/Mkk)^ 
and (Ms/Mkk)^, respectively, expecting that the UV divergence is tamed at least by stringy 
states at an energy scale around Mg. 
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It is too naive, however, to make a guess based only on the superficial degree of divergence. 
Let us take 1-loop correction to a gauge coupling constant as an example. 1-loop correction to 
1/ Qym is expected naively to be of order (M^/Mkk)^, since the superficial degree of divergence is 
+6 for amplitudes with only two gauge fields in the external lines in 9+1 dimensional spacetime. 
The 1-loop correction, however, is known in orbifold calculations and is of order (M^/Mkr)^ 
[27] . This explicit example clearly indicates that the guess based on the superficial degree of 
divergence is too naive. 

1-loop IPI effective action on flat 9+1 dimensional spacetime can explain why the 1-loop 
threshold correction is of order (M^/Mkk)^, not (M^/Mkk)^- If 1-loop IPI effective action had 
a term 1/a''^ tr(|Fp), then its dimensional reduction would give rise to a threshold correction of 
order (R^/a^) = (M^/Mkr)^- But, it is known from worldsheet calculation [28l [29] that there 
is no such correction. Absence of an 0{Mg) term in the 1-loop effective action on fiat 9+1 
dimensional spacetime corresponds to the absence of an 0{Mg) term in the 1-loop threshold 
corrections to the gauge coupling constant of compactified models. On the other hand, the 
1-loop IPI effective action on flat 9+1 dimensional spacetime contains a term [SU] 

d^^y—ti(\F\^)ti(\F\^), (2.91) 
a' 

and dimensional reduction of this term explains the threshold correction of order (M<j/Mkk)^ 
in compactified models. 

String-scale dependence of observables (such as the gauge coupling constant) in compactified 
models may be determined by dimensional reduction of IPI effective action of string theory 
on flat 9+1 dimensions. This prescription is known to work for threshold correction of gauge 
couplings. Although it is desirable to check with orbifold calculations whether this prescription 
is correct for other observable quantities, yet it does not seem to be terribly wrong. Flat 
9+1 dimensional spacetime and its compactification share spacetime structure (including an 
extended supersymmetry) and interactions on it at short distance. Thus, the string scale 
dependence may be the same for both. 

Let us use this prescription for the 1-loop amplitude of the 5-if(5) mass term; a brute 
force world sheet calculation would only be possible for orbifold compactification, and even in 
orbifold, it would take some time. If we are to use known results on IPI effective action on 
9+1 dimensions instead, we have already seen in the context of tree-level contributions that the 
mass term of chiral fermions in question cannot be obtained through dimensional reduction. 
The same is true for the 1-loop IPI effective action as well, since the result at tree-level is based 
only on S0(9, 1) invariance. Therefore, the 1-loop amplitudes with higher Kaluza-Klein states 
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and stringy states in the loop are not likely to contribute to the total amplitude. On the other 
hand, contributions from low-lying Kaluza-Klein states are not captured by the world-sheet 
calculations on a flat 9+1 dimensional spacetime, since particle spectra on compactified space 
are quite different from those on a flat space at that scale. These IR contributions have already 
been estimated earlier in this subsection. Thus, the rough estimate (I2.83P is vahd, providing 
our prescription is justified. 

The same prescription can be applied to the 1-loop contribution to the /i-term in the 4+1 
model that we described in section 12. 2[ Assuming that the non-vanishing gravitino mass is 
from the 3-form flux, we find that no such term can be written down in effective action on 
9+1 dimensional spacetime with a coefficient of positive power of 1/a' = M^. Thus, the 1-loop 
contributions from around the Kaluza-Klein scale dominate the amplitude, and hence it is 
likely that the estimate following from fl2.42p remains valid. 

This prescription is also applied to the R-parity violating coefficients e for the kinetic mixing 
in fl2.74l) . This R-parity violating effect can be erased away from the kinetic terms by field 
redefinition such as fl2.77p and (12.791) . but it reappears in mass matrices and trilinear interactions 
that involve massive states. Certainly no R-parity violating trilinear interactions are generated 
among massless modes of the MSSM after the redefinition. It also turned out in section [273] that 
mixing of massless eigenstates in the external lines are irrelevant to the effective couplings of the 
dimension-5 proton decay operators. This additional R-parity violating mixing, in principle, 
contributes to R-parity violating non-renormalizable operators in the effective theory, once 
the massive states are integrated out. Since we will discuss R-parity violating dimension-5 
operators in section 12.51 we need an order-of-magnitude estimate of the additional R-parity 
violating mixing from the 1-loop amplitudes. UV-finite 1-loop amplitudes in fl2.75l 12.761) take 
account only of contributions where low-lying Kaluza-Klein modes are running in the loop. 
Until now, we have postponed studying whether contributions from higher Kaluza-Klein states 
and stringy states dominate over those from low-lying states. We will now use the prescription 
to make an estimate of the UV contributions. 

The 1-loop kinetic mixing fl2.74p can be obtained through dimensional reduction, if the 
1-loop IPI effective action contains a term of the form 

^ triF,^iD,D,rF!iFp^)d''y (n > 0), (2.92) 

where vector indices 7,a,/3 and 7 are contracted by metric. In fact, there are no terms 
cubic in field strengths in the 1-loop IPI effective action [281 129] . Although the IPI effective 
action has a term quartic in field strengths, the leading 0{l/a') term is totally symmetric in 
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the four field strengths [30], and is factorized into two i^g-singlets, as in (I2.9ip . Since none 
of 5^ N H cannot be factorized into two singlets of U(l)x,g7 C Es, the 1-loop kinetic mixing 
amplitudes are not obtained from dimensional reduction of quartic terms. Thus, there is not 
even a term proportional to a'~^. All other terms with a coefficient in negative power of a' 
do not give rise to the kinetic mixing terms through dimensional reduction. Thus, we expect 
that the contribution from higher Kaluza-Klein states and stringy states is not proportional 
to a positive power of string scale, possibly because of some cancellation due to extended 
supersymmetry of those particles. The R-parity violating coefficients e, therefore, will not be 
much larger than the contributions from low- lying Kaluza-Klein states in (12.751 I^TT^ . and we 
have an estimate 

see comments after fl2.83p for the meaning of (A^). 

R-parity violating vev's deform mass matrices at tree-level as in fl2.36l [27551 12. 58^ r2.60p . and 
the mixing angles are of the order y'^ (N) /Myj^. The additional mixing due to the 1-loop kinetic 
mixing is of order e. Since the estimate (12.930 contains an extra 1-loop factor, this additional 
mixing (I2.77[[2?79l) is quantitatively negligible. Since the 1-loop mixing is qualitatively the same 
as the tree-level mixing, and negligible quantitatively, we will not make a distinction between 
the two different bases of (I2.77[ 12.790 in the following. 



2.5 R-parity Violating Dimension-5 Operators 

Low-energy effective theories are described only by light degrees of freedom, with all heavy states 
integrated out. Numerous non-renormalizable operators are generated as heavy states are inte- 
grated out, and naturally R-parity violation will also manifest itself in the non-renormalizable 
operators. There are four dimension-5 operators that are SU(5)gut invariant and are odd under 
the R parity: 

VTeff. 3 -^10olOolOo^(5)o + -^^(5)ot(5)o^(5)o5o (2.94) 
in the superpotential and 

Keff. 3 -^10,10,% + -^10oi^(5)o^(5); + h.c. (2.95) 

in the Kahler potential^ [2]. We will discuss in this subsection whether these four operators 
are generated. For operators that are generated, the effective mass parameters Meff. of these 



■^■^ Although these operators are trilinear in low-energy degrees of freedom, we only mean renormalizable 
trilinear ones in the superpotential by trilinear R-parity violation. 



37 




Figure 6: Typical Diagrams which generate W 3 101010iJ(5) for the 4+1 model (a) and for 
the 3+2 model (b). 

operators are expressed in terms of parameters of microscopic descriptions such as Mkk, y'^ {N) 
etc. 

2.5.1 Dimension-5 Operators in the Superpotential 

Let us begin with the first operator in (12.941) . This operator is generated by exchanging massive 
states in the SU(5)gut-(5 + 5) representations, as in the super Feynman diagrams in Figure [61 
Since this operator breaks R parity, R-parity violating vev's (A^'^) or (A^) have to be inserted 
either in the internal line or external lines in the diagram. Propagating in the internal line are 
massive states in the SU(5)gut-5 + 5 representations. Following the argument in (12.721) . one 
finds that diagrams with the mixing in the internal line yield 

-^~l/fi)(M-i)5^(5)y[\) (2.96) 

in both the 4+1 and 3+2 model. 

In the 4+1 model, one finds that the contribution from diagrams with the mixing in the 
internal line becomes 

treating as if Mjj were a2x2or3x3 matrix. We have seen in section 12.31 that this finite 
rank treatment is sufficient in making a rough estimate of M^s,, although there are inifinite 
Kaluza-Klein particles. Contributions due to IOq-IO'* in the external lines are of order 

1 ... ... y'i^) <^'> ... ... y'i^) (^^> 
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Figure 7: Examples of super Feynman diagrams which generate W 3 Hq{5)Ho{5)Ho{5)5o in 
the 4+1 modeL 

ff(5)o-5* mixing in the external line yields an amplitude of order ( ]2.97p . If all the massive states 
are around the Kaluza-Klein scale, Mkk, then l/M^s, ~ y'^y'^{y'' (^N^))/M^^. Since l/Meg. ~ 
y^y'^/Miai^ for the dimension-5 proton decay operators (12.631) . it is likely for y" i^N^'^ ^ Mkk 
that the R-parity violating interaction AVT = lOolOolOo-ffo is weaker than the dimension-5 
proton decay operators in the 4+1 model. 

In the 3+2 model, note first of all that mixing in the external lines is irrelevant. Massless 
eigenstates IOq remain pure U(l)5-7-eigenstates, and there is no mixing for these external states. 
The mixing of Hq into U(l)57-eigenstates 5* or does not generate the opeartor /^W = 
lOolOolOo-f/^O; since we have seen in section 12^3] that there is no way to generate either AW = 
lOolOolOoB or AW = lOo IOq IOq i^'' in the 3+2 model. We further notice that there is no 
contribution from mixing in the internal line by using fl2.96l) and a 3 x 3 matrix fl2.73|] . The 
same is true for the case with a > pairs of like states with masses of order y'^ (^)- On 

the other hand, if there are b > pairs of zero-modes H{5)q-5q that become massive in the 
presence of y'^^-^ (A^) 7^ 0, the operator W 3 IOq IOq 10o-f/^(5)o is generated, with 

The second R-parity violating operator in fl2.94p is generated through diagrams in Figure [7] 
in the 4+1 model and Figure [8] in the 3+2 model. There, massive SU(5)GUT-singlets are 
integrated out. 

In the 4+1 model, the effective mass scale from the diagram Figure [7| (a) is approximately 

1 ..^^^yJ^^^JyT^rmiMj^, (2.100) 



Meff. Mg Mr Mkk Mkk Mr 
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where Mr is defined in ( 12.351) . Thus, the operator HqHqHq5o is larger or smaller than the 
other R-parity violating dimension-5 operator lOolOolOo-ffo in the 4+1 model, depending on 
whether Mr < Mkk or not0 if difference among various trilinear couplings and mass spectra 
are ignored. If there are Kaluza-Klein zero modes in (as assumeqf^ in Table [2]), a diagram 
Figure [7] (b) also contributes and 

<i)^(3)-;-^, (2.101) 



which tends to be larger than fl2.100p : see (12.351) . Therefore, in a crude approximation that 
all the Kaluza-Klein towers begin at a common scale Mkk and that y'^ i^^) is somewhat 
smaller than Mkk, R-parity violating AW = HqHqHq^q is enhanced by ((y*^ (A^))/Mkk)~^ 
relatively to the dimension-5 proton decay operators, and the other R-parity violating operator 
AW = lOolOolOo-f^o suppressed by ((y*" (A^))/Mkk)~"^ in the 4+1 model. Here, we assume 
that there are zero modes of A^, and difference among all the trilinear couplings are ignored. 

In the 3+2 model, an amplitude corresponding to the super Feynman diagram in Figure M 
is 



Meff, Mg MR,efr. 

where MR^gfr. is defined in (I2.57p . Thus, we see that an R-parity violating dimension-5 operator 
is generated in the 3+2 model, independent of high-energy spectrum in the SU(5)gut-5 + 5 
sector [i.e., whether 6 = or not]. For cases with 6 > 0, the second term in (12.941) has larger 
coefficient if Mr^^s. < {u^ (^))^/Mkk, and otherwise, the first one has a larger coupling (12.991) . 



2.5.2 Dimension-5 Operators in the Kiihler Potential 

Two R-parity violating dimension-5 operators in the Kahler potential ( 12.95^ are also generated 
through 1-loop diagrams in Figure [9] (in the 4+1 model) and those in Figure [10] (in the 3+2 
model). Only one Feynman diagram is shown for each operator in each model, although there 
are many others. We roughly estimate contribution from a given diagram, by treating Kaluza- 
Klein towers as if they were finite number of massive states in the 3+1 dimensions. Estimates 

■^^If the low-energy neutrino masses come from the see-saw mechanism involving right-handed neutrinos, then 
Mb cannot be larger than lO^^GeV. If A/kk is around A/qut ~ lO^^GcV or even larger, then Mn < AIkk- 

■^^Tliey are necessary if masses of low-energy neutrinos are to be generated from the seesaw mechanism 
involving right-handed neutrinos. Majorana masses from (j2.34p are not sufficient without a zero mode of N 

m- 



40 



(N)(N) 



.5^ 



(AT) 
Ho 



Figure 8: An example of super Feynman diagrams which generate W 3 Ho{5) Ho{5) Ho{5)5o 
in the 3+2 modeL 




Figure 9: One-loop diagrams in the 4+1 model that generate R-parity violating dimension-5 
operators in the Kahler potential. 
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(b) 



Figure 10: One- loop diagrams which generate 10105'^ (a) and 10H{5)H{5y (b) for the 3+2 
model. 



of the UV contributions are postponed until the end of this subsection. This is just like what 
we did in sections [2^ and [2^ for other 1-loop amplitudes. To take the operator K 3 IOqIOoBJ 
as an example, contributions from diagrams in Figure [9] (a) and Figure [10] (a) are roughly of 
order 



4+1 model 



3+2 model 



1 



Meff. 
1 



M 



eff. 



167r2 |M|4 
167r2 |M|2 ' 



(2.103) 
(2.104) 



where |M|'s in the denominators are energy scales where dominant contributions come from 
in integrations of loop momenta in field theory on 3+1 dimensional spacetime. All the con- 
tributions, including those above, are roughly of order 1/Mkk with a 1-loop factor ?/^/(167r2) 
and an extra suppression factor {y" (A^))/Mkk- Here, we assume that all the mass parameters 
appearing in the loop amplitudes are around a common Kaluza-Klein scale Mkk- We do not 
go into detailed discussion of whether those operators have enhanced contributions when states 
with masses of order 0{y'^ i^)) exist. 

The rough estimates above only account for 1-loop contributions from low-lying Kaluza- 
Klein particles. As we discussed in section 12. 4[ however, such contributions may not dominate 
in 1/Meff. in principle. The same operators in effective theories may also have tree-level con- 
tributions. It is also possible that the 1-loop amplitudes are dominated by contributions from 
higher Kaluza-Klein modes and stringy modes. 

First, it is easy to see that super Yang-Mills theory on 9+1 dimensions does not give rise 
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Figure 11: Starting with a trilinear superpotential W 3 AqBqC + y'^CDQ (N), one might think 
that an effective operator K 3 AqBqDI could be generated in the Kahler potential as in the 
diagram (a). However, such a trilinear term does not exist when the state going to the right is 
a massless eigenstate. We are saying that the amplitudes in Figure [7] and [8] with mixing in the 
external lines do not vanish for massless states, while the amplitude in (a) does. To understand 
this subtlety better, it helps to draw Feynman diagrams in terms of mass eigenstates as in (b). 
It is clear that a mass flip is necessary on the external line going to the right in (b), whereas 
it was not necessary in Figure [7] and Figure [HI This means that Kahler terms of the form 
K 3 AqBoDI, with all of external fields being massless, cannot be generated from a diagram 
like (a). 

to the effective operators fl2.95p at tree-level; here, by super Yang-Mills theory, we only mean 
all the interactions that are derived from (Kaluza-Klein decomposition of) the leading order 
action on 9+1 dimensions, {l/{a'^gg)) tr(|Fp)(i^°?/. Supergraphs like those in Figure [TT] (a) exist 
only when one of chiral multiplets in the external lines is a massive mode. If all the external 
states are massless eigenmodes, such amplitudes vanish. For more detailed explanation, see the 
caption of Figure [TT] 

Tree-level contributions to the effective interactions (12.951) may also come from dimensional 
reduction of higher-order terms in a'-expansion in the tree-level effective action on 9+1 dimen- 
sions. Let us take 

4 + 1 model: P 1010 (W) \d = 0^,5'^ dnOlO (N") + ■■■ , (2.105) 
3 + 2 model: P (N)* 1010\d = 8^,5^ (N)* d^lOlO + ■ ■ ■ (2.106) 

in the effective theory on 3+1 dimensions as an example for concreteness; all the fields in the 
right-hand sides of fl2.105[ 12.1061) are complex scalars in given representations. Such interactions 
can be obtained through dimensional reduction, if string theory effective action has a term of 
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the form 



(2.107) 



with indices a-7 and a-7 contracted covariantly by metric haa- It is important that all the four 
field strength tensors are in a single trace, not factorized into two. This is because there is no 
way separating four representations 5~'^, two 10~^'s and 1"^^ in (12.1051) [resp. 5^^, two 10~^'s 
and 1+^ in fl2:T06D ] into two singlets of SU(5)gut x U(1)x C [resp. SU(5)gut x \]{l)q, C ^s]- 
In fact, it is known that the 0{\/a') F^ term of the sphere amplitude is factorized into two Eg 
singlets in the Heterotic Eg x Eg string theory [32]. Thus, the effective operators (12.1051 [2.106P 
on 3+1 dimensions are not obtained through dimensional reduction of tree- level amplitudes, 
at least from 0{l/a') operators. If the effective action on 9+1 dimensions contains a term 
{I/Qs) tT{F^)d^^y, which comes at the next order in a'-expansion, then its dimensional reduction 
may give rise to a largest possible tree-level contribution 



Here, suppression due to small overlap of wavefunctions is ignored, and hence y's and gyu are 
much the same in this crude approximation. Rough estimates of 1-loop amplitudes (12.1031 
12.1041) are 1/Mcs_ ~ (?YM/(167r^) x {gyu i^) /^kk) this level of crude approximation. Thus, 
although tree-level contributions do not have a 1-loop factor (7YM/(167r^), they are suppressed 
by (Mxk/Ms)^ or even more, and it is not very likely that they dominate over the estimates 



Let us now turn to 1-loop amplitudes. Rough estimates (12.103^ 12.1041) take account only of 
contributions with low-lying Kaluza-Klein modes running in the loop. This is why the estimates 
based on field theory on 3+1 dimensions are UV finite; superficial degree of divergence is —4 
and —2 in the 4+1 and 3+2 model, respectively. When all the Kaluza-Klein states run in 
the loop, however, loop momentum is effectively integrated on a 10- dimensional space, and the 
superficial degree of divergence becomes +2 and +4, respectively. Naive guess would be, then, 
to multiply (M^/Mkk)^ and (M^/Mkk)^, respectively, to the rough estimates in (12.1031 CTID . 
because stringy states may set in around the string scale M<,, taming the UV divergence. 

We have already seen in section 12. 4[ however, that this guess is too naive. We introduced 
a prescription in section [2^ for how to guess whether the amplitudes 1/Mes. are enhanced 
by a positive power of (M^/Mkk): read off Ms = l/y/a' dependence of terms in the 1-loop 
IPI effective action on 9+1 dimensions whose dimensional reduction give rise to low-energy 
effective interactions of our interest. This prescription is based on a belief that Calabi-Yau 




(2.108) 



following from (ICT^CTMD . 



44 



compactification and a flat 9+1 dimensional spacetime have a common short-distance spacetime 
and interactions on it, and that compactification introduces only MxK-dependence, but no extra 
Ms-dependence. 

1-loop effective action of the Heterotic Eg x Eg string theory on a flat 9+1 dimensional 
spacetime has been calculated from torus amplitudes. 0{l/a') terms quartic in gauge field 
strength are symmetric for four field strengths for the torus amplitude [30], and are factorized 
into two i?8-singlets, 

-^tr(|F|2)tr d^^. (2.109) 

Dimensional reduction of this term does not give rise to the effective interactions fl2.105[ I2.106P 
for the same reason explained for the tree-level amplitudes. All other operators that appear at 
higher order in the a'-expansion do not have coefficients in negative power of a'. Therefore, 
the 1-loop amplitudes for the effective operators fl2.105[ 12.1061) will not have coefficients in a 
positive power of Mg, and in particular, the estimates (12.1 03^ I2.104p do not need to be changed. 

2.6 Brief Summary 

The minimal supersymmetric standard model has four independent chiral multiplets in the 
(2,-1/2) representation of the SU(2)l x U(l)y gauge group. The dimension-4 proton decay 
problem implies that there must be some structure in the vector space C spanned by the four 
chiral multiplets. The most popular approach has been to introduce an unbroken discrete 
symmetry such as R parity. The vector space splits into two (or more) subspaces that trans- 
form differently under the discrete symmetry transformation. In the case of an unbroken Z2 
symmetry, for example, 

= Span {Hd} © Span {Li, L2, L3} . (2.110) 

Chiral multiplets in the two subspaces have different interactions, and hence down-type Yukawa 
couplings can exist, while trilinear R-parity violating couplings do not. 

An alternative approach presented in this article does not need to assume that a discrete 
symmetry remains unbroken. In order to solve the dimension-4 proton decay problem, the 
vector space £ does not need to split. It is sufficient to assume that C has a subspace 

£d£^; (2.111) 

chiral multiplets that belong to the subspace have restricted interactions. Three independent 
degrees of freedom in £^ are identified with the lepton doublets, and the remaining one degree 
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of freedom in Lj ^ L^"^ with the down-type Higgs doublet. It is not strictly necessary that 
£ has a subspace like Cf^ for H^, that is characterized by its restricted interactions. Suppose 
that a theory has a U(l) gauge symmetry with a negative [resp. positive] Fayet-Iliopoulos 
parameter, and there are chiral multiplets iV'^ with a positive charge and iV with a negative 
charge. Then, N'^ [resp. A^] develops a non-vanishing expectation value, absorbing the Fayet- 
Iliopoulos parameter in the U(l) D-term potential. If there is an interaction (12.341) in the 
superpotential, then the vev of [resp. is set to zero dynamically because of F-term 
potential coming from (12.341) . The U(l) symmetry broken only by positively charged [resp. 
negatively charged] chiral multiplets leaves a structure (12.1111) in £, and similar ones in all the 
other vector spaces of massless modes of given representations of SU(3)c x SU(2)i x U(l)y as 
well. See (12.401 [2391 12.6ip . and discussions that follow. Trilinear R-parity violating operators 
(11.11) are absent, no matter how large the Fayet-Iliopoulos parameter and (A^*^) [resp. (iV)] 
are. At the same time, the interaction (I2.34p gives rise to Majorana masses for right-handed 
neutrinos. Therefore, the absence of trilinear R-parity violating terms of massless modes is 
closely related to the Majorana masses. 

If the Heterotic string theory is compactified with a vector bundle given by an extension, 
then one always has the U(l) gauge symmetry, Fayet-Iliopoulos parameter and interaction 
(I2.34p . Thus, dimension-4 proton decay problem is solved (see footnotes [T3] and [2T1) . and 
Majorana masses are generated for right-haded neutrinos. One does not need to choose a 
vacuum by hand only from special points in the moduli space, since an extra unbroken symmetry 
is not necessary any more. 

Since the U(l) symmetry is broken spontaneously, it does not have a control over the 
Kahler potential of low-energy effective theory. Kinetic mixing terms between Lj and can 
be generated, but they can be erased away by field redefinition, while keeping the structure 
(12.1 lip . See discussion that follows (12.740 for more detail. Thus, the structure (12.1 lip and the 
likes for other representations remain to be a valid solution to the dimension-4 proton decay 
problem. 

The broken U(l) symmetry does not leave a structure like (I2.11ip for massive fields. Mass 
matrices are deformed by the vev's (A^) (either (A^'^) or (A^))- Once massive states are inte- 
grated out to obtain low-energy effective theory, U(l)-symmetry-breaking vev's can enter in the 
denominator of coefficients in the low-energy effective superpotential. Since positively [resp. 
negatively] charged vev's in the denominator supply negative [resp. positive] U(l) charges, 
simple selection rule for the superpotential based on U(l)-charge counting ceased to be valid. 
One cannot claim that an operator with positive [resp. negative] U(l) charges is not allowed 
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in the effective superpotential, even when the U(l) symmetry is brolcen by vev's of positively 
[resp. negatively] charged fields. 

This argument does not mean that there is no rule at all in the low-energy effective super- 
potential. Let us suppose that we have a superpotential of a fundamental theory with a typical 
mass scale M. Find vev's for all the fields, expand fluctuations around the vev's, diagonalize 
mass matrices, and integrate out massive states to obtain low-energy effective superpotential. 
There is nothing special in this process. But, it is important to note that terms trilinear in 
massless fluctuations remain unaffected in the last step of integrating out massive states. Thus, 
the selection rule using U(l)-charge counting can be used to constrain trilinear terms in the 
low-energy effective superpotential. Combined with an observation above that U(l)-breaking 
kinetic mixing in the Kahler potential can be absorbed be redefinition of chiral multiplets while 
keeping the lower triangular nature of mass matrices, we see that the selection rule can be 
used in eliminating the dimension-4 proton decay operators (11. ip . It is not that this idea has 
not been presented anywhere else [SI H El E]- But we consider that this article strengthens 
theoretical basis for using the U(l)-charge selection rule for this purpose, by clarifying the limit 
of the rule as well as some subtleties and logical steps that we tend to overlook. 

When it comes to dimension-5 operators, first of all, there is a good chance that operators 
with negative [resp. positive] U(l) charges are in the low-energy effective superpotential, if 
holomorphic insertion of positively [resp. negatively] charged vev's can make them neutral 
under the U(l) symmetry. Such operators may be in the superpotential of a fundamental 
theory from the beginning, and even if this is not the case, they tend to be generated in the 
process of integrating out heavy particles. 

Dimension-5 operators with n > [resp. n < 0] U(l) charges still have a chance of being 
generated, if some massive states in an appropriate pair of vector-like representations have mass 
parameters carrying m > n [resp. m < n] U(l) charges. Those operators may be generated 
when those states are integrated out. If a dimension-5 operator does not satisfy this condition, 
it is not generated in the low-energy effective superpotential. This observation was used in one 
of the models described in this paper in eliminating the dimension-5 proton decay operators. 

Since R-parity is broken by vev's (A^), it is interesting how R-parity violation appear in 
low-energy effective theory. Particle contents at low-energy and their charges under the broken 
U(l) symmetry alone do not have enough power to answer to this question. This is because, 
for example, we never know whether dimension-5 operators with positive [resp. negative] U(l) 
charges are generated without knowing types of particles around the mass scale M and interac- 
tions that those particles have. Since those operators tend to have couplings enhanced by some 
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positive power of (M/ (A^)), if they are generated, such operators tend to be the most important 
among dimension-5 operators. Thus, we need a (well-motivated) theoretical framework that 
specifies types of massive particles and their interactions. 

We turned to a class of Heterotic string compactification that we mentioned above. Two 
models (called 4+1 model and 3+2 model) that belong to this class were analyzed in this article. 
In the 4+1 model, R-parity violating operator AW = H(5)oH{5) H{5)o5o tends to have an 
enhanced coupling than the dimension-5 proton decay operators. The other R-parity violating 
dimension-5 operator AW = 10olOolOo-??(5)o is also generated with a suppressed coefficient. 
Typical effective couplings of those dimension-5 operators are of order 



where Mkk is the Kaluza-Klein scale and |/ooh's are trilinear couplings among two light states 
and one massive state. Enhancement or suppression factor is roughly given by an appropriate 
power of {y (N) /Mkk)- In the 3+2 model, the dimension-5 proton decay operators are absent 
in the effective theory. On the other hand, both R-parity violating dimension-5 operators can 
be generated^ Thus, the dimension-5 operator in the effective superpotential (apart from 
iJ(5)5iJ(5)5 that we may have already seen) with the largest coupling breaks R parity in both 
models. 

Since the theoretical framework in this section is realized in string theory, we can discuss 
i) R-parity violating operators generated by loop diagrams in the effective Kahler potential, 
as well as ii) those that originate from higher order terms suppressed by gravitational scale 
(or string scale). Bilinear R-parity violation W 3 fiiLiHu comes from a 1-loop correction to 
the Kahler potential, with /ij proportional to the SUSY breaking. We found that tree-level 
contribution is absent when the gravitino mass originates from a 3- form flux. Thus, /ij's are 
suppressed by 1-loop factor and {y (N) /Mkk) relatively to gravitino mass, so that it can be 
very small compared with the electroweak scale. 

Both dimension-5 R-parity violating operators in the effective Kahler potential (12.951) are 
also generated at 1-loop level. At the crudest level of approximation, coefficients of these 
dimension-5 operators are of order 

^ ~ ^ (2 113) 

Meff. IGtt^Mkk' ^ ■ ^ 

where now youns are trilinear couplings among two heavy states and one light state. Each 
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An exception is the operator AW = IOq IOq IOq -ff (5)o for the case with spectrum characterized by 6 = 0. 
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operator is further suppressed or enhanced in an odd power of {y {N)) / My,]^. See section [53] 
for more details. 

Sections 12 . 11 - 12 . 51 are written as if the theoretical framework is based on the Heterotic E^xE'^, 
string theory. But absence of (or negligible) trilinear R-parity violation and presence of all other 
R parity violation rely only on a few ingredients: an extra U(l) symmetry, its charge assignment, 
non- vanishing Fayet-Iliopoulos parameter and interactions fl2.20H2^27|l and fl2.46ti235|) . The 
algebra of Eg alone is sufficient in justifying all these assumptions. Therefore, M-theory dual 
and F-theory dual vacua that share the Eg algebra should also share the same qualitative 
conclusion (c.f. [7]). However, we have also used some properties specific to the Heterotic 
string theory, not just algebra of Eg. That is where we argued that tree- level contributions to 
fl2.74[ 12.801) are absent and those to (12.951) can be negligible compared with those from 1-loop, 
and that the 1-loop contributions are not proportional to a positive power of string scale. Thus, 
one has to study such questions separately for non-Heterotic vacua. 

3 Phenomenology 

Most of phenomenological constraints on R-parity violating couplings that have been discussed 
in the literature are on the trilinear couplings in the superpotential (see [1] for example). Since 
the theoretical framework in the previous section predicts that trilinear R-parity violation is 
either absent or highly suppressed, most of them, including the dimension-4 proton decay prob- 
lem, are no longer a problem. R-parity violating interactions dominantly come from bilinear 
terms and dimension-5 operator in this framework. In the (virtual) absence of trilinear R-parity 
violation, constraints and predictions on nucleon decay are totally different, and some other 
constraints are expressed in much simpler ways because of fewer R-parity violating parameters. 
It is one of the purposes of this section to summarize phenomenological constraints on R-parity 
violating couplings in the (virtual) absence of trilinear R-parity violation. In section 13. 3[ we 
study lower bounds on bilinear and dimension-5 R-parity violating couplings, requiring that 
the LSP in the visible sector decay before the period of BBN. Recent developments in under- 
standing of constraints on hypothetical particles from the BBN are reflected in the analysis. In 
section [231 we discuss nucleon decay processes induced by squark exchange diagrams combining 
dimension-5 and bilinear R-parity violating operators. 

The theoretical framework in the previous section not only predicts which operators are 
generated in low-energy effective theories. It also provided order-of-magnitude estimates of 
those R-parity violating couplings. This means that we can analyze whether or not it can 
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survive various constraints. The order-of-magnitude estimates are compared with a constraint 
from low-energy neutrino masses in section 13. ![ and with those from washout of baryon /lepton 
asymmetry in section [3l2l 

We are by no means the first to study phenomenology of R-parity violation in the absence 
of trilinear R-parity violating couplings in the superpotential. Indeed, there are a plenty of 
literature, especially those focusing on bilinear R-parity violation. Since R-parity violation is 
dominated by bilinear terms at the renormalizable level in low-energy effective theory of our 
framework, technical results in those literature are quite useful. The appendix provides a quick 
summary of such results used in this section. 

3.1 Neutrino Mass 

Bilinear R-parity violation in fl2.86l) introduces neutrino-higgsino mixings, and hence the neutrino- 
neutralino system gives rise to an extra see-saw contribution to the neutrino masses [TU]. As 
the neutrino masses are bounded from above by cosmological observations, there is an upper 
limit on bilinear R-parity violation. Using m,y < leV, the limit is placed on a misalignment 
parameter |e'| [33| : 

\e'\^3.0xl0~' x^xf^y. (3.1) 
cos p V 1 eV / 

|e'| measures the difference between f^i/fio and Vi/vd, its definition is found in [33] (and also 
flA.9l) in the appendix). The limit (13. ip also means that bilinear R-parity violation can provide 
the dominant part of neutrino masses in the case that (13.11) is marginally satisfied. 

Since we have obtained a crude order-of-magnitude estimate of bilinear R-parity violation 
Hi in section [21, it is interesting to see if it satisfies the constraint (13.11) . |e'| in (13. ip is roughly 
of order 0{\ni/v\) in the absence of alignment (see the appendix and discussion below); here, 
V is the electroweak scale. Using (I2.83P in yUj = Cjm3/2, 

In gauge mediation scenario, it is very easy to satisfy the bound if /ij is proportional to gravitino 
mass, but at the same time, it is very unlikely that the neutrino mass comes from bilinear R- 
parity violation. In the case that 7713/2 ~ 100 GeV, the neutrino mass bound (13.11) is also 
satisfied with a safe margin for a reasonable choice y ~ 10~^ and y'^ (N) /My^y. ~ 10~^. In 
this case, however, it is hard to say whether or not the neutrino-higgsino see-saw mechanism 
provides one of dominant contributions to the low-energy neutrino masses, because of large 
theoretical uncertainties associated with (13. 2p . 
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The misalignment parameter | e' | flA.Qp can be smaller than ^/^~JiMJJM)^ for some scenarios 
of SUSY breaking [SU [35]. As we see in the appendix, each misalignment parameter e- (z = 
1,2,3) can be as small as 10"^ x (yUi/yUo) in minimal SUGRA scenario, but it is unlikely that 
they are even smaller than that. Thus, the neutrino mass bound on ^i/v may be relaxed by 
about two orders of magnitude. 

In anomaly mediation scenario, Bi ~ ^3/2 x /ij, so the misalignment parameter |e'| can 
be larger than naive expectation 0{\^i/ ^\). If Bq is somehow of order 0{v'^) and /iq ~ 
then the misalignment parameters e[ are approximately e- ~ —B^/Bq ^ x (fii/v) ~ 

10"''^ X (ni/v). Thus, in this case, the limit on /Xj/f is stronger by three orders of magnitude. 

3.2 Washout of Baryon/Lepton Asymmetry 

Any interactions violating either baryon or lepton number symmetry could wash out the 
baryon/lepton asymmetry of the universe that is once generated after inflation. Thus, such 
couplings should not be too large. All the bilinear R-parity violating operators, 

AW = f^iH^Li (3.3) 

and 

AKoft = -BiHJi + mlQ,H% + h.c, (3.4) 

break lepton number symmetry. When all the dimension-5 R-parity violating operators are 
rewritten in terms of MSSM chiral multiplets, 

[I01010i5^(5)]^ ^ [QQQH,]p + [IJQEH,]^ ^ O3 + O4, (3.5) 

[H{^)hH{h)H{h)\^ ^ -L [H^LH^H,]^ ^ O,, (3.6) 



[10105^ 



eff. L^" '' ^•'-'y ^•'-'iU J D 

1 r^.. 1 



+ + (^10, (3.7) 



Me, [^(5)10^^(5)^^ - ^ - (3-8) 

it is easy to see that each operator breaks either baryon or lepton number sjTiimetry. Requiring 
that they do not wash out baryon/lepton asymmetry, upper bounds on these R-parity violating 
couplings are obtained [36] . 

Tiny trilinear R-parity violation fl2.85[ 12.871) in the A+l model 

- y'it ~ y'it ^^^^^^ ' ^^^^^^ - a + (3.9) 
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also breaks lepton and baryon number symmetry. However, the coupling constants are so small 
for Mh ^ Mgut (^gut ~ lO^^GeV) that they are irrelevant to the washout of baryon/lepton 
asymmetry. 

Bilinear R-parity violating operators are most dangerous at late time, as interaction rates 
scale as oc T while the temperature T is higher than the electroweak scale, whereas the Hubble 
parameter scales as during radiation dominance. Assuming that baryon/lepton asymmetry 
was generated when the temperature was higher than the electroweak scale, and requiring 
baryon/lepton asymmetry not to be washed out, one can obtain upper limits on bilinear R- 
parity violation. Applying the results of [H7] to the case without trilinear R-parity violation, 

mini=i,2,3 ( — ) ^ 10"^ X {y^)'^ ~ 3 x 10"^ x cos/3, (3.10) 

where = y^ is the bottom Yukawa coupling. The rough estimate 03.21) of Hi satisfies this 
constraint easily. Constraints on soft bilinear R-parity violating parameters in [3S] are a little 
stronger: 

(^--\^ 10"^. (3.11) 

As long as Bi ^ 0{v x /Xj), however, the estimate (13.21) also satisfies this constraint from 
washout. 

Dimension-5 operators, on the other hand, are more relevant at higher temperature. Inter- 
action rates of R-parity violating processes scale as F ~ 10~^T'^/Mj^ {i = 3, 4, 6, 7, 8, 9, 10), and 
hence these operators would wash out baryon/lepton asymmetry right after it was generated, 
if they ever would. Requiring that R-parity violating processes caused by the dimension-5 
operators are out of equilibrium at that time, it follows that [36] 

«!=10"G«Vx(^^)*, (3^12) 

where Tab/al stands for the temperature of baryo/lepto-genesisQ Note that (13.121) should be 
satisfied for all seven operators (l3.5H3.8p . If the baryon asymmetry of the universe originates 
from thermal leptogenesis, we know that T^b/al ^ 2 x lO^GeV [39]. In this scenario, all the 
effective mass scales M, in (l3.5H3.8p have to be larger than about lO^^GeV. We have seen in 
section [2] that the zeroth-order approximation of I/M3 4 6 in the superpotential and I/M7 8 9 10 in 
the Kahler potential are y'^/M^K and (?/^/167r^)l/MKK, respectively. Therefore, it is quite easy 

■^^ Limits from the absence of washout of baryon/lepton asymmetry are more comphcatcd if T^s/al ^ 
lO^^GeV, because the sphaleron process is out of equihbrium when T ^ lO^^GeV. 
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to satisfy the constraints from the washout of baryon/lepton asymmetry, if the Kaluza-Klein 
scale is around the GUT scale Mqut- 



3.3 LSP Decay 

Given the order-of-magnitude estimate of bilinear R-parity violation in fl3.2p . it is very unlikely 
that the lightest supersymmetric particle in the visible sector (hereafter we call it the vLSlxj) 
has a lifetime much longer than the age of the universe. If the temperature of the universe was 
once sufficiently high, and thermal relic of vLSP is left after the temperature drops below the 
electroweak scale, then the relic vLSP has to have decayed before the period of BBN. 

Constraints on R-parity violating couplings from BBN were discussed already in [TOl HHl HI] 
and in more detail in other papers that followed. In this subsection, we reanalyze the BBN 
limits in the light of the latest understanding of the impact of hadronic energy injection or 
of stable charged hypothetical particles during the period of BBN. We ignore tiny couplings 
( I2.85P that are present in the 4+1 model, and assume that all the renormalizable interactions 
violating R-parity originate from bilinear terms, as predicted by the framework in section [21 
Limits on dimension-5 R-parity violating operators are also derived. The following study covers 
two typical possibilities: either the vLSP is a bino-like neutralino or a scalar tau (stau). 

The vLSP is no longer a candidate of dark matter. But gravitino with 7713/2 ~ (9(1) GeV 
can be a good candidate of dark matter [32]. Axion and a strongly interacting stable particle 
with a mass of order 100 TeV (if such a particle exists) [43] can also be dark matter. 

Neutralino Let us first begin with the bino-like neutralino vLSP. The thermal relic density 
of the neutralino is very model-dependent. Here we adopt the prediction of the so-called "bulk" 
region of mSUGRA parameter space: 



mpFjo ~ 4 X 10~^° GeV 



100 GeV 



(3.13) 



Nucleons and anti-nucleons produced in jets from the decay of semi-stable particles would 
contribute to p ^ n conversion processes and change the fraction of ^He. Thus, any semi-stable 
hypothetical particle X with baryonic branching fraction Bh has to have short enough lifetime 
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We coin a term vLSP because whether gravitino is hghter than vLSP is sometimes not quite important. 
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Tx ^ 0.1 sec, as long as the relic density satisfies 



{mY)^ xBhZ lO^^'J-lO^'GeV. (3.14) 

If there is enough mass difference between the vLSP neutralino and the LSP gravitino {m^o — 
m3/2 3> Mz), two-body decay x° — ^ + Z is possible at tree level, and ~ Thus 
f l3.14p is typically satisfied in this case, and r^o ^ 0.1 sec would be in conflict with the primordial 
^He abundance. This sets a limit 

m,,.S0^8MeVx|A^,,|x(^)'x(^)% (3.15) 

where N^^ = — sin^vK^i^ + cos is the Z component of the neutralino vLSP. 

If the mass difference is not large enough (m^^o — 7713/2 ^ Mz), however, x° — > ip3/2 + 7 is 
the only two-body decay at tree level, and baryons are not contained in the decay products 
of this dominant decay mode. Baryons are produced only in three-body decay processes, and 
baryonic branching fraction is of order Bh ~ O{10~^). In this case, f l3.14p is not satisfied, and 
hence the limit from p ^ n conversion Tx ^ 0.1 sec does not apply. Deuteron production due 
to hadro dissociation of ^He sets the most stringent limit instead; tx ^ 10^ sec is required for a 
hypothetical particle X if |44] 

{mY)x xBhZ lO^^^GeV, (3.16) 

and this condition is satisfied by the typical relic density of neutralino (13.131) even after multi- 
plying Bh ~ 0{1). Thus, 

™3/.S24MeVx|A^,|x(^)*x(Iti^)*. (3,17) 

where A''i^ = cos 6*1^ A''-^^ -|- smOwNyyyo- 

If gravitino is not as light as specified above, a bino-like vLSP has to decay fast enough 
through R-parity violating operators before the period of BBN. Bilinear R-parity violation 
induces R-parity violating vertices v-x^-Z and e^-x^-W^ in mass-eigenstate basis |16] (see 
also (]A.421{Al44l) in the appendix). If there is enough final state phase space for the two body 



The lower bound in (|3.14p is sensitive to the estimate of the error of measurement of ''He abundance, and 
here we adopt a conservative one in [48]. If the error estimate gets even larger, it brings up the lower bound 
in (|3.14p . but the change is rather slight. On the contrary, if the error estimate gets smaller, the constraint 
Tx ^ 0.1 sec applies to a semi-stable particle X with much smaller {mY)x x Bh- 
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Figure 12: Feynman diagrams for three-body R-parity violating decay of neutralino. Feynman 
rules for the R-parity violating vertices, each shown as a blob, are found in [IHl HZ] and also in 
the appendix of this article, (a) is with a virtual Z/W , (b) with a virtual squark and (c) with a 
virtual non-colored scalar, (c) is only an example among several kinds of similar diagrams with 
a linear combination of 0° = {h^ , , ,ul), 0+ = e|., e^^.) or its complex conjugates 0" 
as the virtual particle. 



decay processes 



V 



Z and x*^ " 
> Z + vjv) 



2 
IGvr 

2 

IGtt 



W"^ {i = e,fi,T), 
m^o gz 



E 

i 



Mz 2 



(3.18) 
(3.19) 



where we have assumed a little hierarchy Mz,w -C Msusy (^susy ~ ^1,2, /^o); for a bino-like 
neutralino vLSP, the ^^.^N^j^ terms in flA.42|-|A.44p are proportional to (M^/Msusy); while all 
other terms in the vertices (]A.42l IA.44 ) arjH to (Mz/Mc^ttsv)^- Using the approximate form 
of ^-.^ in flA.lSP in the appendix and assuming N-^^j^ ~ 1, we find that 



10 



-10.5 



> 



200GeV\2 1 /O.lsec 
X X 

Ml J 10 cos /3 \ 1 

n conversion, because ~ 



pn 



(3.20) 



0{1) in these processes, and 



Here we used the constraint from p 
f l37[3|) satisfies fIXTD . 

If there is not enough final state phase space for x° v + Z^v + Z and l"^ + , there is no 
two-body decay processes at tree level. Examples of Feynman diagrams for three-body decay 
processes are shown in Figure [T2l it is easy to see that the baryonic branching fraction is of order 



''^ Note that there is a partial canceUation between U. ?i~N*~g and 50^1 So)* in (|A.44k see (IA.17I) 

'^d Ifrf '^'^d ^^d 

and (|A.22p . The neutrahno-VF-charged-lepton vertex (|A.43p is negUgible compared with (jA.44[) . because ah the 
coefficients in (|A.43p are suppressed further by mi/yito, where are charged lepton masses; see (|A.23|) . 



55 




Figure 13: Feynman diagrams for R-parity violating radiative decay of neutralino. (a) is an 
example of diagrams with and — 1-charged fermions ip' = {W~ , ,eLk)'^ in the loop. 



(b) is an example of those with charged scalar fields 



and charged fermions ip"^, where 



unity, no matter which one of the diagrams (a)-(c) gives the dominant contribution. Although 
a kinematically allowed two-body decay process x° — > u + '-f does not contain a baryon or even a 
hadron in the decay products, it is generated only at one-loop (Figure [T3|) . and is not expected 
to be much faster than the three-body decay processes. Thus, the baryonic branching fraction 
of R-parity violating decay of neutralino vLSP is of order unity even when m^o ^ mz,w- This 
is in contrast to the case of neutralino decay to gravitino. In order not to spoil the prediction 
of the standard BBN, we neecj^ r^o ^ 0.1 sec, just like in the case of m^o ^ Mz.w- K is likely 
that the diagram Figure [12] (a) gives larger contribution to the decay rate than the others] n 
and 



2 X 3.65 
1927r3 



2x9 
1927r3 



m|o. 



(3.21) 
(3.22) 



where 



g'^/{8M^). Interference between final state neutrinos in fl3.2ip 



is ignored. C^^^ and C^^^ are dimensionless coefficients of Ui-Xi^Z and e^-Xi^W^ vertices in 



To be conservative, the constraint from p ^ n conversion should be replaced by that from excessive 
deuteron production through hadrodissociation, t^o < 10^ sec, relaxing the lower bound on the bilinear R 
parity violation by one order of magnitude and a half. This is because the yield of neutralino in (|3.13p can be 
so small for m^o <, Mz,w that (|3.14p may not be satisfied. It should also be remembered that relic density of 
neutralino vLSP can be even smaller than the estimate (|3.13p in some parameter region of the SUSY breaking. 

The amplitude for Figure [12] (c) is proportional to a Yukawa coupling, and will remain relatively small 
unless tan (3 is very large. 
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( IA.42P and (lA.44p . respectively, and are of the order e[ or fii/fio- Thus, we conclude that 



\/|C(^)|2 + 0.4|C(^)|- 
10-10 



O.lsec 



I pn 



lOOGeV 



(3.23) 



if m^o ^ Mz,w 

The bino-like neutralino vLSP can also decay through dimension-5 operators violating R- 
parity before the period of BBN, even if bilinear R-parity violation is too small to satisfy fl3.20p 
or (13.231) . The decay modes through Og.s are similar to those using bilinear violation, since 
lepton-Higgs mixing is induced when some of Higgs multiplets are replaced by their vev's. In 
the case of m^o ^ M^z, two-body decay processes of a neutralino vLSP u + Z and 

— > £^ + W^. are induced by Oq. The decay widths for these modes can be calculated by using 
Goldstone equivalence theorem, that is, treating the longitudinal component of gauge boson in 
the final state as Goldstone boson, with mass Mw or Mz and coupling equal to that of Higgs 
boson. These two decay modes are comparable, and the decay width is 



r(x° 



Z) 



IGvr 



On the other hand, the only two-body decay process induced by Og is x° 
mixing is induced only for charged leptons. The decay width is 



(3.24) 

9± + ]yT since 



c;0 



r(x 



IGvr 



A^^^o cos /3 



(3.25) 



Applying the constraint r^^o ^0.1 sec, we have 

Me < lO^^^GeV x (10 cos /5) 



-1/2 



lOOGeV/ 



Ms < 5 X lO'^GeV x (10 cos /3) 



lOOGev) 



1/2 



(3.26) 
(3.27) 



where we assumed that m^o 



and 



Mb 



sin 6w sin /5, 



Mz_ 

Ms 



sin 9w cos /?. 



(3.28) 



In the case of m^o ^ M^^z, a vLSP neutralino decays to three fermions with a virtual gauge 
boson in the intermediate state. Oq with two H^s replaced by their vev's induces a mixing in 
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Figure 14: Examples of Feynman diagrams for neutralino decay that involve dimension-5 R- 
parity violating interactions, (a) uses O4, and (b) Og. 

the neutral part H^L'^, and consequently the R-parity violating couplings x'^-u-Z and 
W^. We thus find that the constraint (13.231) leads to 

The dimension-5 operator contributes to three-body decay of vLSP neutralino only through 
virtual-VT processes. We find, after a calculation similar to what we have had so far, that the 
upper bound of Mg is 

M, , :0"OeV X COCO.,) {?^) (^)' i^f. (3.30) 

Feynman diagrams in Figure [TH show that there are three-body decay processes for neu- 
tralino vLSP that involve dimension-5 R-parity violation 4 in (13.51) or 07,9.10 in (13.71) . The 
decay width of these processes are given by (c./. [49] ) 

r(x° - ff* - fff") - M\N,s\ \aofo J 2n T' (3-31) 

192(27r)^cos''t^VK "^0 

where Yf is the hypercharge of a fermion /, and ttlq the mass of a virtual sfermion /. Effective 
coupling Aeflf is given by (if^) /M3 4 for C3 4 and by m^/My 9 10 for 07,9,10, where mj- is a mass 
of an outgoing fermion (either /' or /", depending on which is the heavier one). Combinatoric 
factors such as the number of colors and final states are ignored here because it depends on 
which operator is dominant and what kind of flavor structures they have. Since a pair of 
quark and anti-quark is always in the final states in these three-body decay processes, baryonic 
branching fraction is of order unity. Imposing the limit r^o ^0.1 sec as before, we find that 

o 1 /100GeV\V "^0 \2/0.1sec\^ 

Aeff > 2.3 X 10-«-— — — ^ . 3.32 

\N^^\Yf V m^o J VlTeV/ V ^pn J 
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This limit corresponds to 

M3, ^ lO'GeV X (lOcos,) (^)' (3.33) 



^7,9,10 

Only either one of (13.331) and (13.341) has to be satisfied. 

To summarize, vLSP neutralino must decay fast enough not to spoil BBN. For m^o S> Mz 
[m^o <^ Mz], neutralino decays to gravitino fast enough for gravitino mass satisfying (I3.15P 
(13.171) ]. For larger gravitino mass, at least either one of R-parity violating couplings must be 
large enough; either one of (13.201) [ (I3.23p ]. (I3.33P and (13.341) must be satisfied. 

Stau Next we consider the case where the vLSP is stau. Typical thermal relic of stau is 

2 



100^) ' ^^-^^^ 



.lOOGeV/ 

and if it decays too late, successful predictions of the BBN are no longer valid. 

Problems with the BBN can be avoided for sufficiently light gravitino, as the two-body 
decay t ^ r + ?/'3/2 can be fast enough for stau to decay before the period of BBN. The most 
stringent constraint on late-time decay r — > ?/;3/2 + r is to require [50l [5T] that the primordial 
value of the abundance ratio {n6u/n7u)p remains unchanged by the stau-catalyzed process 
in the presence of long-lived stau [52]. The upper bound is roughly ^ 10^ secj^ which 
corresponds to m3/2 ^0.1 GeV for ~ 100 GeV. 

For gravitino mass m3/2 ^ 0.1 GeV, stau must decay through R-parity violating operators 
instead. This requirement sets a lower bound on the R-parity violating couplings. A stau 
decays dominantly to a pair of leptons in the presence of bilinear R-parity violation. Although 
it also decays to a pair of quark and anti-quark, such an amplitude is proportional to tau lepton 
mass, and the branching fraction is suppressed by a factor of order {rrir tan /^/ttisusy)^, which is 
not expected to be larger than about 10~^. See [53] or the appendix for more. Since the typical 
yield of thermal relic of stau in (I3.35P with Bh ^ 10~^ does not satisfy (13.161) . the deuteron 
constraint r ^ 10^ sec does not have to be imposed for stau vLSP decay. The most stringent 
limit on the stau lifetime comes from the stau-catalyzed process, and we need to require that 



^■^ The stau-catalyzed process is a problem in the range m^rlV 10 ^^GeV. Typical thermal relic of stau 
vLSP (|3:35l) is within this range. 
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< 10^ sec. Using flX45U07|) Fl 



for k = 1,2, 



(3.36) 



and we find that 



\ 



k=l,2 ^'^^ 



in-i3/in /'lOOGeVWlO^sec 
^ 1.8 X 10 ^''(10 cos P) [ — 1 



(3.37) 



Stau may also decay fast enough to a quark-ant i- quark pair through dimension-5 operators 
O4 10 with 



qq) 



167r 



(3.38) 



even if bihnear R-parity violation is not large enough to satisfy fl3.37p . The effective couplings 
Acflf for this decay are Acg. ~ {H^) /M^^ for and ~ nif/Mio for Oiq. If the two-body decay 
through O4 or Oio is to be the solution to the BBN problem, the effective coupling of these 
operators has to be large enough, so that r^^ ^ 10^ sec; since the baryonic branching fraction 
of the R-parity violating decay through O4 or Oio is of order unity, fl3.16p is satisfied. We thus 
find that 

lOOGeVy /I02sec^ 



Aeff. > 3.3 X 10"^^ 



^d~had 



(3.39) 



This lower bound on Acg. is equivalent to either 

M4 ^ 5 X lO^^GeV X (10 cos /3) x 



'^d—had 



100 GeV/ V 102 sec 



or 



Mio ^ 3 X 10'^ xnif X 



12 



/ 



lOOGeV/ VlO^sec 



Td~had \ 2 



(3.40) 



(3.41) 



It is sufficient if either one of these is satisfied. 



tion 



he dimension-5 operator Cg induces stau decay through effective bilinear R-parity viola- 
^ On the condition that ^ Vd, dominant contribution comes from W 3 (vJ^/Mq) L^H^ 



Decay amplitudes are dominated by diagrams involving the vertices (jA.46p . as long as we assume large 
tan f3. But the vertex ()A.45|) can be more important for small tan f3. 

Although Oq also provides effective trilinear R-parity violating vertices with only one of Higgs fields replaced 
by its vev, their contribution to stau-vLSP decay is negligible in the limit where left-right mixing is ignored, 
since we assume that the vLSP r is the scalar partner of r^. 
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which causes lepton-Higgs mixing on 
mixing gives dominant contributioijl^ 



y in the neutral part. A decay mode using bino-neutrino 



r(r 



1 



16 



sin P sin 9w 



(3.42) 



and the constraint from BBN is given by 



Mfi < 4 X 10 



13 



-V 



T"6Li? 



(3.43) 



lOOGeV/ V200GeVy V200GeVy VlO^sec. 
The dimension-5 operator also contributes to stau decay through effective bilinear R- 
parity violation, with one of Higgs fields replaced by its vev. In contrast to the case with Og, 
lepton-Higgs mixing occurs only in the charged part. The decay width of stau vLSP depends 
on whether gravitino mass is large or not. In addition, if 3> Mz + M-\y, a two-body decay 
process r ^ — Z + is also possible, without any of Higgs fields replaced by their vev'sF^ 

To summarize, if the vLSP is stau, and if 7713/2 ^ lOOMeV, one of the R-parity violating 
couplings should be sufficiently large, so that stau can decay fast enough and the standard 
predictions of the BBN are not spoiled. It is sufficient if either one of R-parity violating 
operators allow stau-vLSP to decay. 

The neutrino mass bound (13.11) and the absence of washout of baryon/lepton asymmetry 
(13. lot 13.121) require that none of R-parity violating couplings are too large. On the other 
hand, some of R-parity violating couplings should be sufficiently large so that the vLSP can 
decay before the period of the BBN, unless gravitino mass is sufficiently small. Figure[T7| (a)-(d) 
summarize the relation among these phenomenological constraints in simplified two-dimensional 
parameter space characterizing the order of magnitude of bilinear and dimension-5 R parity 
violation. 

Flavor dependence of R-parity violating couplings has been ignored so far. Let the effective 
operator O4 in (13. 5p has a flavor dependent coefficient W 3 C\%QiUjEkHd/ M^^, for example. It 



is E. 



1-^(4)12 



that determines the decay rate of stau vLSP, while what matters to the washout 



Another decay mode using down-type Higgsino-neutrino mixing can give comparable contribution only if 
tan/3 is very large. 

"^^ It is interesting to note that the branching fraction of various decay modes stau vLSP vary so much, 
depending on which R-parity violating operator is responsible primarily for the vLSP decay. Oq predicts that 
stau-vLSP decays dominantly into final states including r, whereas the branching fractions of decay modes to r 
and /U or e have fixed ratio if the conventional bilinear R-parity violation dominates |53j . A branching fraction 
to /X or e can be larger that that to r, if Og dominates. If R-parity violating decay of stau vLSP is observed at 
the LHC, branching fractions of various decay modes provide valuable information on physics behind R-parity 
violation. 
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of lepton asymmetry is the smallest among {k = 1, 2, 3); baryon/lepton asymmetry 

is washed out completely, only when all three B/3 — {k = 1,2,3) symmetries are broken 
by some interactions in the thermal equilibrium. Therefore, with flavor structure of R-parity 
violation, the true allowed region in the parameter space tends to be wider than it appears in 
Figure [171 

3.4 Nucleon Decay 

Dimension-5 R-parity violating operators in ( 13.51) and Oj in (13.71) break baryon number 
symmetry, and hence proton may decay. In the 4+1 model, baryon number is broken also 
in the trilinear R-parity violation Oq in f l2.85p . In section [3.4.11 we derive limits on R-parity 
violating couplings from proton lifetime. We further discuss in section 13.4.21 how to probe 
R-parity violation through nucleon decay experiments. 

3.4.1 Limits on R-Parity Violating Couplings from Proton Lifetime 

Decay products of a proton contain at least one fermion, and it must be one of e, /i, u or their 
anti-particles if it is a particle in the Standard Model. Thus, proton decay is induced by squark 
exchange diagrams combining two R-parity violating operators; one breaks baryon number and 
the other lepton number. O3, Oj and Oq are candidates for the former, while renormalizable 
interactions u^jj^ePL/Rd in (]A.36p . d*^^j^e'^PL/jiu in (1A.37I) and q*uq + h.c. in (lA.38f[Al4T]) for 
the latter. The interactions (lA.36f[X.4ip originate from bilinear R-parity violation. R-parity 
violating dimension-5 operators Odin fl3.5p and Og^io in fl3.7p can also play the same role for 
necessary lepton number violatioiO- 

Feynman diagrams in Figure [15] show squark-exchange diagrams combining baryon-number 
violating O3, Or and Oq and a lepton-number violating (]A.37p . Diagrams using (1A.37P con- 
tribute to nucleon decay processes such as p ^ 7r° + e"*" where a positively charged lepton is in 
the decay products. After integrating out SUSY particles, effective operators become 

Dimension-5 operators Oe.s also violate R-parity and lepton number. They induce effective bilinear R- 
parity violating operators when some of Higgs fields are replaced with its vev. The consequence can be discussed 
by borrowing the constraint on bilinear R-parity violation (|3.54p . but it turns out that such constraint is looser 
if Mq^s are assumed to be in the same order as M4_9_io, and hence we do not discuss them here. 
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Figure 15: Feynman diagrams leading to proton decay p M" + £"*" (and neutron decay 
n — >■ M~ + that originate from the combination of two operators Qq- llA.STp in (a), 
f [07l) in (b) and O'^- I^A^ in (c). 
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Os-(^ £off. = {aU.^w-) i^LdL){uLeL), (3.44) 

Or-i^KM C,^. = :^Jf {aU^LW-) dMR^^UL){uLeL): (3.45) 

dL 

Generation indices are dropped, and the CKM matrix elements that appear in (1A.37P are 
ignored, dj^j^ [resp.ex/^] in this subsection just means down- type quarks [resp. charged leptons], 
not just the down quark [resp. electron] in the first generation. U ^^^^^ and U~^g- are matrix 
elements appearing in diagonalization of charged fermion mass matrix, and are proportional 
to bilinear R-parity violation. For their definition, see [IHl HTj or the appendix of this article. 
Hadronic matrix elements 

{M\p}\uL{uLdLMk)) = W{q^)uL{k), (3.47) 
{M\p)\uL{uRdRMk)) = W\q^)uL{k) (3.48) 

have been obtained by lattice simulation [SU ES] for neutral pseudoscalars M° such as vr^ and 
and the form factors W and W are of order Aq^^j^, where Aqcd ~ 300 Me V is the QCD 
scale. Although the matrix element 

(M'^muLd^idRa^'uLMk)) = W"{q^)uL{k), (3.49) 

is not known to our knowledge, we would expect that the form factor W" is of order Aq(jj3. 

We have discussed in the previous section how the effective mass scales of dimension-5 
operators depend on parameters of microscopic descriptions such as the Kaluza-Klein scale 
Mkk and a suppression (sometimes enhancement) factor associated with a U(l) symmetry 
breaking {y {N) /Mkk)- Writing down the three contributions to the proton decay amplitudes 
in terms of those parameters, we find that they are proportional to 

—iau^w) - y'^^^igu.^ ), (3.50) 

Mg^^ e^M/ ^ y Mkk Mkk 
m/^^^-^-^ 167r2^MKK Mkk^^^^^"^' ^ ^ ^ 

y'^J^iy'^U ~ ) ^ jiLj-^'^iy'^U ~ ) f3 52) 



where we ignored a possibility that some of the dimension-5 operators are enhanced when there 
are vector-hke pair of particles whose masses are of order y (N) . Rough ratio between these 
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three operators io 



OOH; 



/ COS /3 \ ^ 




10-^ 


V 0.1 J 







X (2/1 



OHHj 



/ cos/3\ 




10-2-5 


V 0.1 )\ 







X (?/ooh)(2/ohh)(2/hhh) 



0.1 

cos P ^ 
(3.53) 



where we used expressions for U~^y^^ and U^^g- in (lA.22p . and assumed that y'^ ~ fng/vd (m^ 
is the s quark mass). Thus the combination Oii- l\A.37\i is hkely to contribute the most. It 
should be remembered, however, that the order-of-magnitude estimates in section [2] have very 
large uncertainties. 

Experimental limits on proton lifetime are roughly around r(p + £"^) ^ lO'^^-lO'^^yrs. 

for major decay modes. This limit can be used to set an upper bound on the bilinear-dimension- 
5 R-parity violation. Using the partial amplitude from the combination 0^- (\A.37h . we finc@ 



l(PGeV\ 
Ms ) 



10 



-9 



2 

m~ 

dL 



;iTeV)^ 



10GeV\ I gV2Mw/M2 cos(3^ 



Vd 



10 



-1.5 



(3.54) 



A similar constraint is obtained for the partial amplitude from Oy- (\A.37\\ when M3 is replaced 
by M7 and Vd by Aqcd- Unless different partial amplitudes cancel one another, each of these 
constraints have to be satisfied, and the one with O3 gives stronger constraint if M3 and M7 is 
in the similar order. 

For even smaller bilinear R-parity violation, proton decay amplitudes are dominated by 
squark-exchange diagrams combining two R-parity violating dimension-5 operators. Typical 
Feynman diagramed are found in Figure [TBI (a)-(c). Requiring that each partial amplitude is 



Sec (|2.112l [2.113P for the meaning of variatfons of Yukawa coupling constants j/ooh, yoHH and ?/hhh. 

For consistent comparison with constraints from the vLSP decay and others, we evaluate effective mass 
scales like M3 renormalized at weak scale here. 

The 1-loop amplitude of Figure [TBI (a) has a larger contribution than a tree-level one with a propagating 
Higgs boson replaced by its vev. The 1-loop diagram is logarithmically divergent, and the divergence is cut off 
at the energy scale of masses of heavy particles that are already integrated out in section [^751 1-loop numerical 
factor including logarithm and (l/167r^) is about 0.4, thanks to the large logarithm. 
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small enough, we find that 

10^^GeV \ /l(FGeV\ 
M3 ) V Af4 ) 



(3.55) 

10^^GeV \ / 10^^GeV \ /10GeV\ / 10-"-^GeV \ 

10^5GeV\ /10i5GeV\ ^ /10GeV\ /10-°-^GeV\ lO'^ 



'10 



) ^ A (^) (^) r ^ (3.-) 



where Oifi is a mixing angle between left-handed and right-handed squarks. 

If gravitino is lighter than a proton, then gravitino can be the fermion in the decay products 
of a proton. A squark exchange diagram in Figure [16] (d) induces proton decay, by combining 
baryon number violating O-s and 

Oa = ^-^D^^^,a^rd,G, (3.58) 

where G is the Goldstino field and (g, ipg) is a pair of complex scalar and chiral fermion of any 
one of chiral multiplets of the MSSM. The effective operator corresponding to the diagram is 

0;-0^ C,s. = /7T \. {d,GLa''a^dL)d,{uLdL). (3.59) 

Relevant hadronic matrix element is of the form [56] 

{M+{p)\a''a>'dMuLdLMk)) = pWuM, (3.60) 
and the form factor W" will be of the order of Aq^^j^. Thus, we should require that 

or otherwise the partial amplitude from O^-O^ would predict proton decay faster than the 
experimental bound. 

Upper bounds on dimension-5 and bilinear R-parity violation derived from proton decay 
are shown in Figure [TTJ along with constraints from neutrino mass (section 13. ip . washout 
of baryon/lepton asymmetry (section 13.21) and BBN (section l3.3p . The four parameter-space 
plots in Figure [T7| are intended just to provide a big picture of the allowed parameter space 
of bilinear-dimension-5 R-parity violation. As the effective mass scale M^s. of dimension-5 
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(a) bino vLSP, m^o 3> mz,w, ■"^■3/2 ^ IMeV, 

(b) bino vLSP, m^o ^ mz, ?^3/2 ^ 30MeV, 

(c) stau vLSP, m3/2 ^ 100 MeV, 

(d) sufficiently small m3/2 

Figure 17: Simplified picture of allowed parameter space of R-parity violation. |e'| parametrizes 
bilinear R-parity violation, and Meg. ~ 4 6,7,8,9,10 sets the scale of dimension-5 R-parity 
violating operators. Hatched areas are excluded. We assumed that T^b/al ~ 10^° GeV and 
vLSP's are thermally produced. 
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R-parity violation, we adopt M3 4 of the operators (^3^4 for demonstration. If M3 4 5,7,8,9,10 
are assumed to be in the same order of magnitude, the operators (^3,4 are the most effective 
pair for proton decay. For washout constraint, there are no significant dependence on which 
operators are used. On the other hand, some of dimension-5 R-parity violating operators among 

7,8,9,10 let the vLSP decay faster than through ^3,4. In such cases, the constraint from BBN 
on dimension-5 R-parity violation is looser than shown in Figure [T71 but general perspective of 
allowed parameter region is not largely altered. 

If baryon asymmetry was generated when the temperature is below the electroweak scale, 
then the upper bound on the bilinear R-parity violation is replaced by the neutrino mass bound 
from cosmology. The allowed parameter space is extended to larger |e'| by about two orders of 
magnitude. 

Flavor structures are not taken into account in Figure [T71 For the case of stau vLSP, allowed 
parameter space may extend to arbitrary small |e'| with Mcs, ~ lO^^GeV, if flavor structures 
of dimension-5 operators are taken into account; see the comments at the end of section 1X51 

The theoretical framework in section [2] provides rough estimates of bilinear and dimension- 
5 R-parity violation, |e'| ~ O{10~^) and Meg. ~ (9(Mkk)- Here Mkk ~ ^gut in a scenario 
where SU(5)gut symmetry breaking is associated with the compactification of spacetime. It 
roughly corresponds to the upper-right corner of the parameter space that survives all the 
phenomenological constraints discussed in sections 13.1113.41 in the big-picture parameter space 
in Figure [T71 

3.4.2 Probing R-Parity Violation with Nucleon Decay 

No evidence for proton decay has been found so far. However, once nucleon decay is discovered 
by experiments, then data of branching fractions of various decay modes can be used to study 
physics behind the nucleon decay. 

Nucleon decay has been discussed in the literature mainly in the context of unified theories. 
GUT gauge boson exchange and colored Higgsino exchange predict dimension-6 and -5 operators 
in the Kahler potential and superpotential: 



All these operators happen to preserve B — L. Therefore, B — L number is preserved in all the 
nucleon decay processes p[n\ M^[M~] + and p[n] — > M+[M°] + u which are predicted by 
conventional (SUSY) GUT's. It is also known [57] that the B — L number is preserved in all 



K 3 E^QU^Q + D^LU^Q, 
W 3 QQQL + UUED. 



(3.62) 
(3.63) 
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Figure 18: A Feynman diagram ioi B — L violating neutron decay. In this example the combi- 
nation of O'^- llA.SGI) is used. 



the baryon-number violating dimension-6 operators of the Standard Model that preserve the 
gauge group SU(3)c x SU(2)i x U(l)y. 

A Feynman diagram in Figure [18] shows that a, B — L breaking nucleon decay n — > M"*" + i~ 
is possible in the presence of R-parity violation. Squark-exchange diagrams involving u*ed 
vertex (1A.36P induce neutron decay with a negatively charged lepton in the final state. This 
B — L breaking neutron decay, coming from a partial amplitude 0-?- ()A.36p . is rather a robust 



prediction of R-parity violation. As one can see in Table HJ all but one combination of two 
operators for nucleon decay allow the B — L breaking neutron decay in the third row|f§ Even 
in gauge mediation scenario with very light gravitino of 7723/2 ~ (9(1-10 eV), the amplitudes 
involving is at most comparable with other decay modes, so n — -|- i~ may be observed 
with a considerable fraction. 

The B — L violating nucleon decay is possible because all the effective operators that we 
obtain after integrating out squarks are dimension-7 or higher. Either Higgs boson (vev) or 
derivatives are involved in those operators. For example. 



Os-m^ C,^.^\y-^^^^d^{lq){qH,), (3.64) 

Tfl— IVl'i 



is a dimension-7 operator in the Standard-Model fields, with the coefficient l/m~ coming from 
squark exchange. After a Higgs doublet is replaced by its vev, this operator effectively becomes 
a 4-fermion dimension-6 operator with a coefficient suppressed only by a single power of a large 



Although proton decay processes p + ly also break the B — L symmetry, there is no way to confirm 

in experiments whether the missing particle is v or something else. 
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03-(ren.) 


07-(ren.) 




03,7-09 


O3-O10 


O7-O10 




p 


n 


M^[M-]+i+ 










(LR) 






p 


n 


M+[MO] + /0 
















M+ + £- 






(LR) 


(LR) 


(LR) 


(LR) 





Table 4: Check marks in this table show that a combination of operators in a given column 
contributes to nucleon decay processes specified in the three rows. M^'" are mesons with given 
electric charges, £^ = e^,/!"*^, and /° a neutral fermion such as v^v and G. "(ren.)" in the 
first two columns mean renormalizable interactions listed in ( ]A.36MA.41[) . Some nucleon decay 
partial amplitudes are generated with an insertion of left-right mixing in the virtual squark 
propagator (see Figure fT6] (c)). Such cases are indicated in the table by "(LR)". 
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Figure 19: Feynman diagrams for K 3 5^ 5^ 5^ 10, (a) for the 4+1 model and (b) for the 3+2 
model. 

energy scale M3. Thus, such an effective operator can be more important than the nucleon 
decay operators 03.62^ 13.631) of conventional unified theories. 

The B — L violating neutron decay is regarded as a signal different from prediction of 
conventional unified theories, and certainly R-parity violation is one of possible explanations. 
One should keep in mind, however, that an effective dimension-6 operator 

K 3 5^5^PlO = L^D^D^Q + D^D^D^E + L^L^D^U (3.65) 

also breaks B — L symmetry, while R-parity is conserved. It is true that this operator cannot 
be generated by gauge-boson exchange processes in any unified theories on 3+1 dimensions, 
as such processes would result in two anti-chiral multiplets and two chiral multiplets. But, 
it will be possible to come up with a model of SUSY GUT with some new matter multiplets 
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Figure 20: A Feynman diagram ioi B — L violating neutron decay through K 3 5^ 5^ 5^ 10. 

and their interactions in the superpotential, so that this operator is generated@ The first two 
operators of the MSSM in fl3.65p generate dimension-7 effective operators in the standard model, 
after sfermions and gauginos are integrated out (Figure [20]) . Therefore, just an observation of 
B — L violating neutron decay is not enough to conclude that R-parity is not preserved. It 
is an interesting question whether it is possible to distinguish nucleon decay due to R-parity 
violation from one through fl3.65p . using data on branching fraction of various decay modes. 
But this question is beyond the scope of this paper. 

It should be possible to derive predictions on flavor dependent branching fractions by ex- 
ploiting flavor structure of R-parity violating operators. All three chiral multiplets Q in 
cannot be in the same generation; one chiral multiplet Q in the second generation is necessarily 
involved. On the other hand, operators flA.36l - [A.41|) and Oq does not introduce a genera- 
tion mixing (except for predictable small flavor mixing from CKM matrix elements in (1A.361 
IA.37P ). Thus, dominant decay modes will include a K meson in the final state, if the nucleon 
decay is dominated by partial amplitudes (9';i- flA.36tiA.4ip or 0^,-0 Quantitative predictions, 
however, are not covered in this article. 

3.5 Brief Summary 

In section [3l we studied phenomenological constraints on R-parity violation, specialized to 
cases when there are only bilinear and dimension-5 R-parity violation. The (virtual) absence 
of trilinear R-parity violation is a prediction of the theoretical framework in section [21 BBN 
constraints were reanalyzed in section 13.31 where we exploited the latest understanding of the 

This operator exists even in the 4+1 model and 3+2 model that we explained in section [2?2l See Figure [T9l 
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impact of hadronic energy injection and presence of stable charged particle during the period 
of BBN. Section 13.41 placed a limit on bilinear-dimension-5 R-parity violation through proton 
decay. An allowed parameter space of bilinear-dimension-5 R parity violation is presented in 
Figure [TTl where constraints from BBN and proton decay are shown along with those from 
cosmological neutrino mass bound and absence of washout of baryon/lepton asymmetry at 
high temperature. Theoretical framework in section [2] predicts a theoretically likely region in 
the parameter space (with very large uncertainties). The region is roughly around the upper 
right corner of the space that is still allowed by all the phenomenological constraints discussed 
above. 

Because the likely region is around the upper right corner, the framework in the section [2] 
prefer stronger R-parity violation within the allowed parameter region. This means that there 
is a chance that R-parity violation is confirmed by experiments. R-parity violating decay of 
the LSP in the visible sector (vLSP) may be observed inside the detectors of the LHC for 
large bilinear R-parity violation. Studies of this signal of R-parity violation are found in the 
literature, and we have nothing to add in this article. The R-parity violating decay of vLSP 
may not be observed in accelerators, however, for smaller bilinear R-parity violation and/or 
light gravitino in the gauge mediation scenario. 

As we discussed in section 13. 4[ nucleon decay can be an alternative way to probe R-parity 
violation. Especially, non-vanishing branching fraction oi B — L breaking neutron decay n — ^ 
M"*" + £~ is a robust prediction of bilinear-dimension-5 R-parity violation, n —>■ M"*" + i~ is 
always predicted except in the decay mode into gravitino, and the decay ratio into gravitino 
can be at most comparable with other decay mode, even if gravitino mass is very small as 
(9(1-10 eV). This is a notable feature of bilinear-dimension-5 R-parity violation because the 
observation of n — M~^+£~ enables us to distinguish it from conventional (SUSY) GUT's, which 
predict only B — L preserving nucleon decay processes. Therefore, nucleon decay experiments 
can be complementary to the R-parity violating decay of the vLSP in the detectors of the LHC 
in probing R-parity violation. 
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A Notes on Bilinear R-parity Violation 



In the MSSM with bihnear R-parity violation, superpotential is given by 

W = yt^U.Q.H^ - yt^DkQjHd - yl^LuE^H^ + /ioi^n^d + ^^^HuL,, (A.l) 
and soft SUSY breaking potential by 

Koft = rnl^ijQlqj + rn]jijdfd''^ + ml^.uf^^^ (A.2) 

+ mloo\hd? + + (A.3) 

+ BoKhd + BihJi (A.4) 

+ A'^y'^u'^qK- AYd^qhd- A^yHe^hd. (A.5) 

In bilinear R-parity violating scenario, there exists a basis of chiral multiplets (if^, i^i,2,3) such 
that the R-parity violation appears only in bilinear terms as above. This basis is called bilinear 
basis. We basically follow the convention of 

We will assume that all the R-parity violating parameters, /ij, Bi and rnj^Q^, are small 
compared with /xq, Bq and mj^QQ and respectively. Anything that are quadratic in the 

R-parity violating parameters are ignored in this article. Sneutrino vev's in the bilinear basis 
are given approximately by [31] 

^ + (A.6) 

mh, + ^ cos(2/?) 

where Vi = (z/j) = (L°), = {H^), f„ = and tan/? = Vu/vd, so Vi < Vd. We do 

not distinguish \/v'^ + vf from Vd, as the difference between them is quadratic in the small 
R-parity violating parameters. Using the leading-order part of the minimization condition 

("^ioo + /^o) + ^cos(2/5) =5otan/3, (A.7) 

the expression above can be rewritten in a useful form |35j : 

Vi Bi - cot /? (/io/Uj + m| ) 



Vd Bo - cot p{nl + m|oo - m| ■ 
Misalignment parameters are defined by \33^ 



(A.9) 

/iO Vd ' 



^^The misalignment parameter ^ in [33|, [T] corresponds to |e'|, norm of e[. Many references introduced 
different notations for these misalignment parameters. To name a few, Ref. |35| introduces ai = Vi/vd, 
7j = /Ji/no^ Si = Bi/Bo, and difference between arbitrary two out of a, 7 and 5 is basis independent. For 
example, (7j — ai) = e-. Reference [46] uses Aj = e-/xofd. C iii [1] is norm of {5i — ai). 
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and its cot (3 expansion is given by 

< = ^ f ) (l + cot 4) - f cot - + cot 

+ ofl^coep^^coep). (A. 10) 

Vy"o -Bo / 

In minimal SUGRA mediation scenario, the first three terms in (lA.lOp vanish in the UV ini- 
tial condition, and the remaining terms are suppressed by tan~^ (3. Thus, e[ can be much 
smaller than Hi/no or Bi/Bq. For large tan/3, however, the first two terms become larger when 
renormalized at the electroweak scale, as they are generated through radiative corrections in- 
volving bottom Yukawa couplings [Ml [35] : there is still a cancellation of about 10~^tan^/3 in 
(yUj/yUo — Bi/Bo) and about 10~^tan^/3 in (?ti|qq — Tn\^^). In the end, can be as small as 
10~^ X 0{iii/ Ho) ^ 10"^ X C(i?i/i?o), but it is unlikely that e-'s are even smaller than that. 

A.l Mass Matrices 

Neutralino— Neutrino Mixing Because of f j 7^ and /ij 7^ 0, mass matrices of neutralinos 
and neutrinos are mixed up. In the gauge-eigenstate basis ip^ = {B, W^, H^, H^, I'i)'^, the 
neutralino-neutrino mass matrix C 3 —%lj^'^M^,^^^ip^/2 becomes 

V "''N;3x4 '"ss / 

where rriss is the contribution from the see-saw mechanism involving right-handed neutrinos, 
and 

/ Ml -g'vd/V2 g'vjV2 \ 

M2 gvd/V2 -gvjV2 

-g'vd/V^ gvd/V2 -/io 

V g'vu/V2 -gvu/V^ -^o / 

{-g'vi/V2, gvi/^/2, 0, -/i, ) • (A.13) 



(A.12) 



7V;3x4 



Gauge-eigenstate basis ip^ and mass-eigenstate basis ip^ are related by ip^ = Nj-^j ■ ip^, where a 
unitary matrix iVyxy makes A^yxr ' ^n-7x7 ' ^7x7 diagonal. Four mass eigenvalues are supposed 
to be around the electroweak scale or SUSY scale, and three others are very small. We will 
write these mass eigenstates as 2,3,4 ^'i,2,3- Ignoring russ and anything that comes at the 
second order of bilinear R-parity violation, the diagonalization matrix is expressed in a form 
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where ^ ■ M^^.^^^ = ?tI;v-3x4 ^^^^ order, and N^^^Mj;;j.^^^N^^^ is diagonal. The suffixes of 
^Q,/3 run over a G {i^i, z^2, ^'s} and /3 G {5, if^, We will keep the lower-right 3x3 

part of A^7x7 as Isxs, so that the three neutrino-like "massless" eigenstates 2^1,2,3 become i^e,^t,r 
approximately. 

Assuming a little hierarchy between the electroweak symmetry breaking scale and the SUSY- 
breaking scale, i.e., Mz^w ^ ^1,25/^05 simple expressions for ^ai3 are obtained. 

Mz 

iu^B - ^sin^vycos/?e^, (A.15) 



Mz 

^o^w^ - -^cos^M/cos/Je^, (A.16) 



/i, + motan/3e^ /a 
e..HO - , (A.17) 

- -^4 (A.18) 

where 

(Mzcos/5)2(Micos2^iy + M2Sin2^H/) 

"^0 = T-T^TT • (A. 19) 

M1M2 

The tree-level contribution to the low-energy neutrino mass from the higgsino-neutrino mixing 
is given by moe'^e'^ [33] . 

Chargino— Charged-Lepton Mixing Charginos and charged leptons also have a mixed 
mass matrix. In the gauge-eigenstate basis ifj^ = (W~^ , H:^ , e'^)'^ and = {W~ , ,eLi)'^ , 
the mass matrix Af^.g^^ in £ 3 —ip'"^ Mg.^^^ip~^ + h.c. is given by 

/ M2 gv^ 

^C;5X5 = 9Vd /iO -y^V, \ . (A.20) 

\ fii y-Vd 

Mass eigenstates = {xi,2^ ^1,2,3) ^~ = (Xi,2 5 1,2,3) are unitary transforms of the original 
gauge eigenstates; 

^- = f/.^-, {i)+ = v-i)^, (A.21) 

where unitary matrices U and V make U*{MgM~^)U^ and V{M~Mg)V~^ diagonal. 

Under the same approximation as above, eigenvectors for the {e^, /i^, Tij-like mass eigen- 
states cl 1,2,3 are 

U V2Mwcos(3 . 
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Eigenvectors for the {e^, fi^, r'^}-like mass eigenstates 2 3 



_ V2Mw{M2 sin p + fip cos p) m'r , 



m 



(A.23) 



Charged Scalar Mixing There are eight complex scalar fields with +1 electric charge in 
the MSSM, namely, if^*, e^, •, one of which is the longitudinal component of W^. It is 
convenient to adopt the Vi = basis, {H'^, L'^). The Vi = basis and the bilinear basis are 
related by 



L' 



1 Vi/Vd 
-Vi/Vd IsxS 



Hd 
L, 



In the new basis (j)^ = (i^"*", H '*, e^, e'^J"^, mass matrix in £ 9 — 0+'l"Mc;^ 



C;8x8 



/ 










\ 




B'o 


Ml 


m^fi'^tanP 


ML, 








rripfip tan P 


Ml,, 


mixpj 




V 




Mho 




Mh, 


1 



where 



5^,^fio + ^M^sin(2/3), 

B[ ^ Bi - {v,/vd)Bo, 

/i • ~ /ij - {vi/vd)iJ,o = ei/Lto 



(A.24) 
is given by 

(A.25) 



(A.26) 
(A.27) 
(A.28) 



(i.e., {B'-, are {-Bj,/ij} in the Vi = basis), and 

1 



Mf. 



Ml 



Ml, 

Ml,, 
ML, 



m 



Mz cos (2/3) + cos^ /3 + /io + yu: 



LOO + ^^1 cos(2/3) + M^, sin^ /? + 



12 

i 1 



Ml 



mixig 



^ m'4. + (^-Mi - J cos(2/?)5,,- + + (m,^)^5. 
= m|p^+(M|-M^)cos(2/3)5p, + (m;)V 

= - Atotan/5). 



(A.29) 
(A.30) 

(A.31) 

(A.32) 
(A.33) 
(A.34) 



Here m'Lo^ '^'liQi ^'Lj ^^"^ ^'lij non-holomorphic SUSY-breaking mass-square parameters 
in the f , = basis. The eigenvector corresponding to the would-be Goldstone boson is oc 
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-Vd, 0, 0, 0, 0, 0, 0). After removing this mode, the mass matrix becomes 

(A.35) 







m^/ig/ COS P 


5,7 COS /3 




nipfip/ cos/5 


Ml,, 


mixpj 




V 5.7 cos /3 




MU 



In the hmit that charged lepton masses rrf are ignored, both M^;^^^ and ?Ti^/i^ vanish, and 
hence e*'s do not mix with the charged Higgs boson and left-handed charged sleptons. 

A. 2 Three Point Vertices 

We present some of three point interactions of the MSSM with bihnear R-parity violation. 
Three point interactions involving one SUSY-particle-like mass eigenstate and two Standard- 
Model-particle like mass eigenstates are relevant to two-body decay of the visible-sector LSP, 
and to nucleon decay amplitudes. 

Squark Yukawa Couplings Here is a list of all the three point couplings that involve a 
squark and two Standard-Model fermions. They are found in [17] . These interactions are 
combined with C3 (13.51) . Oj fl3.7p . Oq (13.91) to generate nucleon decay amplitudes. 
The three point interactions involving charged leptons are 



+ iype,,H-ryS''''^UekPRd,) + {y^V,.HiV:^'''')uUekPLd,) + h.c, (A.36) 

HyW^lHiy^i''"'rdU^'^PRU,) + (yl^t/,,,5-)V;?™^^.(^fci'LU,) +h.c.(A.37) 

Here, ULi and V^^^^dij are complex scalar fields in the chiral multiplets {un, V^^^dLj), and 
ur and dR [resp. = u'^ and d*^ = d'^] in the anti-chiral multiplets W and [resp. in the 
chiral multiplets U and D]. Four component notations ePid and ePud mean e'^di and CLd^, 
respectively, and c'^Plu = u'^Plc and c'^Pru = u'^Prc correspond to clUl and e'^u'^, respectively. 
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Those involving neutrinos are 

A£ = -V2(^^^,^^ + ^^^^^o^uU,.,PLU.)-iyr^^^^^^^ (A.38) 

{j^o.B - f e.,iyo) dUukPLd^) - {yUo,H-J*dUykPRd^) + h.c, (A.39) 
-V2 [-\9'io,^ uU^kPRU.) - {vKo^HofiU^kPLU^) + h.c, (A.40) 
-V2 [+\9'i,^^ dU^kPRdd - {ytio,H-)dU^kPLdi) + h.c. (A.41) 
The four component spinor = {^^ki'i-'k) is Majorana. 

Neutralino Three-Point Vertices Here, we hst three-point vertices involving the neutralino- 
like mass eigenstate that allow its two-body decay to two Standard-Model fields. 

A£ = --^ua^Z.i^ X ^ [^,^~^N,^,.,^ + i,^y,oN + 2i,^^,N .^...^.^ + h.c(A.42) 

-T%-a^W;t\ X -| [-V2V,.^^,Nl~, + l-,.|,.Ar;^„) + h.c, (A.43) 

In the last line, a runs over {5, H^", if^, if^}. Higgs-neutralino-neutrino three point vertices 
are omitted. 

Stau Three Point Vertices When the left-right mixing is ignored, is a mass-eigenstate 
T itself. In this limit, R-parity violating vertices including stau, which can be used for the 
calculation of stau-vLSP decay, are given as 

AC = {V2g'C,jT%UkPRr) (A; = e,/i,r) (A.45) 

+ m,,Hi}^\^kPLT) + {ylUl^^^)r{u,PLe,) {k + r) (A.46) 

-f(l/3^(f/;^. + i,^^.^)r\vrPLr) + h.c. (A.47) 
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